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TIECHNICAL  NARRATIVE 


SUMMARY 


During  FY  92,  work  was  performed  in  five  project  areas  summarized  briefly  as  follows.  A  more 
detailed  account  is  provided  in  subsequent  sections. 

(a)  The  Role  of  Bubbles  in  Ocean  Acoustics.  Significant  research  progress  was  accomplished 
in  the  latest  project  period.  This  progress  is  documented  in  the  format  proscribed  for  rcpori*  t<j 
the  OcccTi  Acoustics  Program  at  ONR,  and  is  appended.  One  detailed  publication  is  also 
attached. 

(b)  Acoustically  Active  Surfaces.  This  research  began  in  1988  with  an  effort  to  develop  smart 
acoustically  active  coatings  using  piezorubber  and/or  PVDF.  ITie  work  with  smart  coatings  has 
expanded  to  include  three  aspects.  First  concerns  the  development  of  actuators  and  sensors, 
especially  those  for  use  at  low  frequencies;  the  second  concerns  the  study  of  different  control 
algorithms;  and  the  third  involves  the  study  of  the  coupling  between  the  active  surface  and  the 
medium.  The  work  in  these  three  areas  is  discussed  in  this  report. 

(c)  Propagation  Physics.  This  research  addresses  the  physics  of  stochastic  and  deterministic 
sound  propagation  and  scattering  in  the  ocean.  The  principal  investigator  was  not  at  the  National 
Center  for  Physical  Acoustics  after  January  1,  1992.  However,  during  the  first  six  m.onths  of 
1992  the  continuous-wave  part  of  the  research  was  completed  and  the  results  were  written  and 
submitted  for  publication  in  the  Journal  of  the  Acoustical  Society  of  America.  A  preprint  of  the 
article  serves  here  as  the  final  report,  as  previously  agreed  in  the  proposal. 

(d)  Transducer  Development.  Most  of  the  year  was  spent  on  transducer  development  and  in 
completing  the  Ph.D.  Dissertation  of  Dehua  Huang;  however,  some  of  the  time  was  spent  in 
completing  research  supported  in  previous  years  under  Navy  Programs.  These  are  included  in  this 
final  report. 

(e)  Graduate  Fellowships.  The  fellowship  program  was  developed  with  the  hope  that 
outstanding  undergraduates  would  be  identified  and  attracted  to  the  University  of  Mississippi  for 
specialization  in  acoustics  at  the  National  Center  for  Physical  Acoustics.  We  believe  that  this 
program  has  given  more  visibility  to  acoustics  as  a  specialization  in  physics  and  engineering,  and 
that  visibihty  is  in  the  best  interests  of  the  Navy.  The  report  on  the  success  of  this  program  as 
well  as  the  work  of  the  five  students  supported  during  FY  92  follows. 


The  Role  of  Bubbles  in  Ocean  Acoustics 


The  Role  of  Bubble  Clouds  in  Low  Frequency  Scattering  from  the  Ocean  Surface 
Principal  Investigator:  I^awrence  A.  Crum 

Research  Accomplished  in  1992: 

Significant  research  progress  was  accomplished  in  the  latest  project  period.  This  progress 


is  documented  in  the  format  proscribed  for  reports  to  the  Ocean  Acoustics  Program  at  ONR,  and 
is  appended  to  this  communication.  One  detailed  publication  is  also  attached. 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Papers  Published  in  Refereed  Journals 


/.  Articles 

1.  “The  underwater  sounds  of  precipitation,”  L.A.  Crum,  R.A.  Roy  and  A.  Prosperetti, 
Naval  Research  Reviews,  XLIV,  2  (1992). 

2.  'Sonoluminescence  and  bubble  dynamics  for  a  single,  stable,  cavitation  bubble,"  D.F. 
Gaitan,  L.A.  Oum,  C.C.  Church  and  R.A.  Roy,  J.  Acoust.  Soc.  Am.  91(6),  3166-3183, 
1992. 

3.  "Low  frequency  scattering  from  submerged  bubble  clouds R.  Roy,  W.  Carey,  M. 
Nicholas,  J.  Schindall  and  L.  A.  Crum,  J.  Acoust  Soc.  Am.  2993  (1992). 

4.  "Acoustically  forced  oscillations  of  air  bubbles  in  water:  Experimental  results",  R.  G. 
Holt  and  L.  A.  Crum,  J.  Acoust.  Soc.  Am., 91, 1924  (1992). 

5.  "Collective  oscillations  of  a  bubble  cloud  as  a  source  of  underwater  ambient  noise  in 
the  ocean",  W.  S.  Yoon  and  L.  A.  Crum,  J.  Acoust.  Soc.  Korea,  10, 47  (1991). 


II.  Proceedings 

1.  "Acoustic  backscattering  from  a  cylindrical  bubble  cloud  in  water",  S.  W.  Yoon,  K.  J. 
Park,  L.  A.  Crum  and  R.  A.  Roy,  Proceedings:  Western  Pacific  Acoustic  Conference.  FV, 
462  0991). 

2.  "Low  frequency  scattering  from  resonant  bubble  clouds",  R.  Roy,  W.  Carey,  M. 
Nicholas,  and  L.  Crum,  Proceedings  of  14th  ICA.  2,  C3-5,  Beijing,  China,  August,  (1992). 


3.  "Underwater  acoustic  backscatter  from  a  cylindrical  bubble  cloud",  S.  W.  Yoon,  K.  J. 
Park  and  L.  A.  Crum,  Proceedings  of  the  14th  ICA.  2.  C3-1, Beijing,  China,  August, 
0992). 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Papers  in  Refereed  Publications 
(cont.) 


4.  "The  acoustic  emissions  from  bubble  plumes  generated  by  an  impacting  water  jet",  A. 
Kolaini,  C.  C.  Church,  R.  A.  Roy  and  L.  A.  Crum,  Proceedings  of  the  14th  ICA.  2,  C3-4, 
Beijing,  China,  August,  0992) 

5.  "Enhancement  of  hydrodynamic  flow  noise  radiation  by  the  regulation  of  air  bubbles 
in  a  turbulent  water  jet",  M.  S.  Korman,  R.  A.  Roy  and  L.  A.  Crum,  Proceedings  of  the 
14th  ICA.  1,  B6-l,Beijing,  China,  August,  0992). 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Papers  Submitted  to  Refereed  Journals 

I.  Articles 

1 .  "Low  frequency  underwater  sound  generation  by  an  impacting  transient  cylindrical 
water  jet",  A.  Kolaini,  R.  Roy,  L.  Crum  and  Mao  Yi,  J.  Acoust.  Soc.  Am.,  (submitted  for 
publication). 

2.  "Bubble  production  by  capillary-gravity  waves",  A.  Kolaini,  et  al.,  J.  Acoust.  Soc. 
Am.  (submitted  for  publication). 

3.  "The  production  of  high-frequency  ambient  noise  by  capillary  waves,"  R.  A.  Roy,  A. 
R.  Kolaini  and  L.  A.  Crum,  In  Natural  Physical  Sources  of  Underwater  Sound,  ed.  by  B. 
Kerman,  Kluwer  Acad.  Pub.,  Dordrecht,  Netherlands,  in  press. 

4.  "Collective  oscillations  in  a  bubble  cloud,"  S.  W.  Yoon,  K.  J.  Park,  L.  A.  Crum,  M. 
Nicholas,  R.  A.  Roy,  A.  Prosperetti  and  N.  Q.  Lu,  In  Natural  Physical  Sources  of 
Underwater  Sound,  ed.  by  B.  Kerman,  Kluwer  Acad.  Pub.,  Dordrecht,  Netherlands,  in 
press. 

5.  "Sound  emissions  by  a  laboratory  bubble  cloud",  M.  Nicholas,  R.  Roy,  L.  Crum,  H. 
Oguz  and  A.  Prosperetti,  J.  Acoust.  Soc.  Am.,  (submitted  for  publication). 


11.  Proceedings 

1.  "Sound  scattering  from  microbubble  distributions  near  the  sea  surface,"  W.  Carey  and 
R.  A.  Roy,  Proceedings:  The  SACLANT  Ocean  Reverberation  Symposium.  Kluwer 
Acad.,  Dordrecht,  Netherlands,  in  press. 

2.  "Low  Frequency  resonance  backscatter  from  near-surface  bubble  clouds",  L.  A.  Crum, 
J.  A.  Schindall,  R.  A.  Roy,  and  W.  M.  Carev,  Proceedings  of  the  SACLANT  Cpttfecen^ 
on  Low  Frequency  Active  Sonar  (to  be  published). 


Lawrence  A.  Crum 
Ronald  A.  Roy 

Books  or  Chapters 


none 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Technical  Reports 

1.  "Lake  Seneca  Bubble  Scattering  Quick  Look",  W.  Carey  and  R.  A.  Roy,  NUWC  Tech. 
Report  (in  press). 


Lawrence  A.  Crum 
Ronald  A.  Roy 

Invited  Papers 


1.  "Low-frequency  backscatter  from  dense  submerged  bubble  clouds,"  R.A.  Roy,  W. 
Carey,  M.  Nicholas,  L.A.  Crum  and  J.  Schindall,  presented  at  the  special  session  on 
Acoustical  Studies  of  Upper  Ocean  Processes  at  the  123’’^^  meeting  of  the  Acoustical 
Society  of  America,  Salt  Lake  City,  UT,  May,  1992. 

2.  "Laboratory  experiments  on  bubble  cloud  oscillations,"  L.A.  Crum,  R.A.  Roy,  M. 
Nicholas  and  A.  Kolaini,  presented  at  the  special  session  on  Acoustical  Studies  of  Upper 
Ocean  Processes  at  the  123*^  meeting  of  the  Acoustical  Society  of  America,  Salt  Lake 
aty,  UT,  May  1992. 

3.  "Surface  scattering  from  microbubble  distributions  near  the  sear  surface,”  W.  M. 
Carey,  R.  R.  Goodman  and  R.  A.  Roy,  presented  at  the  SACLANT  Ocean  Reverberation 
Symposium,  La  Spezia,  Italy,  May  1992. 

4.  "The  acoustic  emissions  from  bubble  plumes  generated  by  an  impacting  transient 
water  jet,"  A.  R.  Kolaini,  C.  C.  Church,  R.  A.  Roy  and  L.  A.  Crum,  presented  at  the 
meeting  of  the  14^^  ICA,  Beijing,  China,  September  1992. 

5.  "Low  frequency  scattering  from  resonant  bubble  clouds",  R.  Roy,  W.  Carey,  M. 
Nicholas,  and  L.  Crum,  presented  at  the  meeting  of  the  14ih  ICA,  Beijing,  China, 
August,  (1992). 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Coritributcd  Papers 

1.  "Air  bubble  interaction  with  a  submerged  axisymmetric  water  jet,"  S.  M.  Cordry,  R. 
A.  Roy  and  L.  A.  Crum,  presented  at  the  122*^^  meeting  of  the  Acoustical  Society  of 
America,  Houston,  TX,  November,  1991. 

2.  "Sound  velocity  measurements  in  a  bubble  cloud,"  S.  A.  Cheyne  and  R.  A.  Roy, 
presented  at  the  122^^  meeting  of  the  Acoustical  Society  of  America,  Houston,  TX, 
November,  1991. 

3.  "Collective  oscillations  in  fresh  and  salt  water  bubble  clouds,"  R.  A.  Roy,  M. 
Nicholas,  K.  Markiewicz  and  L.  A.  Crum,  presented  at  the  122^^  meeting  of  the 
Acoustical  Society  of  America,  Houston,  TX,  November,  1991. 

4.  "Low-frequency  acoustic  scattering  from  a  submerged  bubble  cloud,"  R.  A.  Roy,  L.  A. 
Crum,  M.  Nicholas,  J.  Schindall,  W.  A.  Carey,  W.  A.  Konrad,  W.  J.  Marshall,  E.  C. 
Monahan  and  A.  Prosperetti,  presented  at  the  122^*^  meeting  of  the  Acoustical  Society  of 
America,  Houston,  TX,  November,  1991. 

5.  "Collective  oscillations  of  a  bubble  cloud  as  a  source  of  underwater  ambient  noise  in 
the  ocean,"  R.  A.  Roy,  S.  W.  Yoon,  M.  Nicholas,  K.  J.  Park  and  L.  A.  Crum,  presented  at 
Intemoise  91 .  Sydney,  Australia,  November,  1991. 

6.  "Acoustic  backscattering  from  a  cylindrical  bubble  cloud  in  water,"  S.  W.  Yoon,  K.  J. 
Park,  L.  A.  Crum,  and  R.  A.  Roy,  presented  at  the  Western  Pacific  Acoustic  Conference 
IV.  Sydney,  Australia,  November,  1991. 

7.  "A  laboratory  experiment  on  nonlinear  coupling  between  surface  and  volume  modes 
of  a  bubble,"  Y.  Mao,  L.  A.  Crum  and  R.  A.  Roy,  presented  at  the  123‘^<^  meeting  of  the 
Acoustical  Society  of  America,  Salt  Lake  City,  UT,  May  1992. 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Contributed  Papers 
(cont) 

8-  "Enhancement  of  hydrodynamic  flow  noise  radiation  by  the  regulation  of  air  bubbles 
in  a  turbulent  flow,"  M.  S.  Korman,  R.  A.  Roy  and  L.  A.  Crum,  presented  at  the  123*^^ 
meeting  of  the  Acoustical  Society  of  America,  Salt  Lake  City,  UT,  May  1992. 

9.  "Acoustic  emissions  from  toroidal  bubbles",  A.  Kolaini,  M.  Nicholas  and  L.  .A.  Crum, 
presented  at  the  123rd  meeting  of  the  Acoustical  Society  of  America,  Salt  Lake  City,  UT, 
May,  1992. 

10.  "Acoustic  characteristics  of  laboratory  breaking  waves",  A.  R.  Kolaini  and  L.  A. 
Crum,  presented  at  the  124th  meeting  of  the  Acoustical  Society  of  America,  New 
Orleans,  LA,  October,  1992. 

1 1.  "Sound  scattering  from  near- surface  bubble  clouds,"  J.  Schindall,  W.  M.  Carey,  R.  A. 
Roy,  L.  A.  Crum  and  M.  Nicholas,  presented  at  the  124^h  meeting  of  the  Acoustical 
Society  of  America,  New  Orleans,  LA,  October  1992. 

12.  "Acoustic  phase  velocity  measurements  in  a  bubbly  liquid  using  a  fiber  optic  laser 
Doppler  velocimeter,"  S.A.  Cheyne,  C.T.  Stebbings  and  R.A.  Roy,  presented  at  the  124^ 
meeting  of  the  Acoustical  Society  of  America,  New  Orleans,  LA,  October  1992. 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Conferences  Attended 

1.  Acoustical  Society  of  America,  Houston,  TX,  4-8  November,  1991  (LAC,  RAR) 

2.  Acoustical  Society  of  America,  Salt  Lake  City,  UT,  11-15  May,  1992  (LAC,  RAR) 

3.  ONR  SRP  Status  Meeting,  Victoria,  BC,  19-21  May,  1992  (LAC) 

4.  SACLANT  Symposium  on  Sea  Surface  Sound,  La  Spezia,  Italy,  25-29  May,  1992 
(LAC,  RAR) 

5.  International  Conference  on  Acoustics,  Beijing,  China,  7-15  September,  1992  (LAC) 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Patents  Awarded 


Patent  Applications 


Lawrence  A.  Crum 
Ronald  A.  Roy 


none 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Honors  and  Awards 


•  Ronald  A.  Roy  was  elected  to  Fellowship  in  the  Acoustical  Society 
of  America. 

•  Lawrence  A.  Crum  was  awarded  a  Significant  Alumni  Achievement 
Award  by  Ohio  University. 


Graduate  Students 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Name:  Christopher  Hobbs 

Citizenship:  USA 

Date  of  Graduation:  August,  1992 

Type  of  Thesis:  Principally  Experimental 

Thesis  Title:  Propagation  of  Sound  through  a  Bubble  Screen 

Thesis  Objective:  To  determine  the  acoustical  properties  of  bubble-filled  geometrical 
shapes:  Bubble  screens  are  used  in  certain  applications  to  "screen"  the  noise  produced  by 
naval  vessels  from  radiating,  by  introducing  a  barrier  that  is  nearly  impenetrable  to  sound 
propagation.  We  determined  that  the  screen  itself  may  resonant,  however,  through  the 
collective  oscillations  of  the  bubbles  contained  within  the  screen. 


Name:  Kenneth  Maridewicz 
Citizenship:  USA 
Date  of  Graduation:  August,  1992 
Type  of  Thesis:  Principally  experimental 

Thesis  Title:  Collective  Oscillations  of  Bubble  Columns  and  Screens 

Thesis  Objective:  To  determine  the  collective-oscillation  frequencies  of  specified 
geometries  of  bubble  columns  and  screens:  Clouds  of  bubbles  can  resonate  at 
frequencies  characteristic  of  the  geometrical  shape  and  acoustic  properties  of  the  cloud, 
which  reduce  the  natural  frequencies  to  values  much  below  that  of  the  individual  bubble 
resonances.  Mr.  Schindall  will  study  this  problem  and  its  application  to  surface  scattering. 


Name:  Jeffrey  Schindall 
Citizenship:  USA 

Estimated  Date  of  Graduation:  May,  1995 
Type  of  Thesis:  Both  experimental  and  theoretical 
Thesis  Title:  Scattering  from  Bubble  Clouds 

Thesis  Objective:  To  determine  the  scattering  characteristics  of  bubble  clouds;  It  has 
been  demonstrated  that  bubble  clouds  can  resonate  at  frequencies  that  are  associated  with 
their  collective  modes.  These  frequencies  can  be  quite  low  and  depend  upon  the  total 
volume  of  air  in  the  cloud  and  the  geometrical  shape.  Thus,  these  objects  can  act  as  false 
target  or  be  employed  as  decoys  or  for  calibration  of  active  sonars. 


Graduate  Students 
(cont.) 


Lawrence  A.  Crum 
Ronald  A.  Roy 


Name:  YiMao 
Citizenship:  PRC 

Estimated  Date  of  Graduation:  August,  1993 
Type  of  Thesis:  Both  theoretical  and  experimental 
Thesis  Title:  Free  oscillations  of  Gas  Bubbles  in  Liquids 

Thesis  Objective:  To  determine  the  acoustical  properties  of  gas  bubbles  at  reduced 
pressures  and  elevated  temperatures:  It  has  been  proposed  that  the  surface  waves  on 
bubbles  could  nonlinearly  couple  into  volume  oscillations,  thus  converting  a  nonradi?iing 
energy  source  into  a  radiating  one.  This  research  is  investigating  the  conditions  under 
which  this  process  is  upheld. 


Name:  Sean  Cordry 
Citizenship:  USA 

Estimated  Date  of  Graduation:  May,  1994 

Type  of  Thesis:  Both  theoretical  and  experimental 

Thesis  Title:  Noise  production  by  Bubble  Fission  and  Fusion. 

Thesis  Objective:  To  determine  the  conditions  under  which  bubbles  coalesce  and  break 
up,  and  the  noise  emissions  associated  with  this  behavior.  When  bubbles  are  created  in 
the  ocean,  they  emit  sound,  and  tiius  contribute  to  the  ambient  background;  it  is  thought 
that  adult  bubbles  are  relatively  silent.  However,  if  the  total  energy  of  a  bubble  is 
changed  by  bubble  fission  or  fusion,  then  they  can  radiate  this  additional  energy.  The 
conditions  under  which  adult  bubbles  become  noise  are  being  examined. 


Lawrence  A.  Crum 
Ronald  A.  Roy 

Postgraduate  Students 

Name:  Michael  Nicholas 
Citizenship:  USA 

Status:  Completed  project;  now  employed  at  NRL 

Work  Objective:  To  undertake  a  systematic  study  of  collective  oscillations  of  bubble 
clouds:  The  concept  of  collective  oscillations  of  bubble  clouds  was  not  demonstrated 
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Transitions 

1.  Dr.  William  Carey,  DARPA,  (703)  696-2339:  Dr.  Carey  has  expressed  an  interest  in 
applying  the  results  of  our  studies  of  bubble  cloud  scattering  to  more  applied  projects 
involving  the  development  of  decoys  and  calibrated  targets.  An  SBIR  has  been  issu^  by 
DARPA  and  we  are  collaborating  with  Quest  Integrated,  Inc.,  to  develop  these  decoys 
and  targets. 

2.  Mr.  Tom  Warfield,  ONT,  (703)  696-5121;  Mr.  Warfield,  a  Program  Manager  at  ONT, 
is  planning  to  fund  us  to  work  with  NRL  to  extend  our  scattering  m,easurements  of  bubble 
clouds  to  a  more  general  understanding  of  low  frequency  acoustic  scattering  from  the  sea 
surface.  We  will  be  working  with  Dr.  Fred  Erskine  (202)  767-3149  and  others  in  the 
NRL  group  to  undertake  these  experiments  in  Exuma  Sound  later  in  1993. 

3  Dr  Sev'ik,  DTRC,  (410)  227-1335;  Dr.  Sevik,  an  Associate  Director  at  David 

Taylor  has  expressed  a  keen  interest  in  our  studies  of  transmission  through  bubble 
screens.  This  facility  pioneered  the  concept  of  bubble  maskers  and  are  interested  in  our 
discoveries  of  the  low  frequency  emissions  that  might  occur  for  these  screens— due  to 
collective  oscillations.  We  are  preparing  a  proposal  for  their  consideration. 

4.  Dr.  Michael  Nicholas,  NRL,  (202)  767-3149;  Dr.  Nicholas  was  a  Postdoctoral 
Research  Associate  under  our  direction  and  is  now  employed  in  the  Underwater 
Acoustics  Branch  of  NRL.  He  will  be  assisting  in  the  analysis  of  the  CST  data  and  will 
likely  transition  insight  gained  in  our  laboratory  from  basic  research  studies  relating  to 
bubble  cloud  scattering  directly  to  his  more  applied  work  at  NRL. 

5.  Dr.  Paul  Hwang,  Quest  Integrated,  Inc.,  (206)  872-9500;  Dr.  Hwang  collaborated  with 
us  in  the  preparation  of  a  SBIR  proposal  to  DARPA  in  which  calibrated  acoustical  targets 
and  decoys  are  to  be  develop^  based  on  the  theory  of  the  collective  oscillations  of 
bubble  clouds.  Later  phases  of  this  work,  if  funded,  will  involve  the  construction  of  these 
targets  for  use  in  the  fleet.  This  particular  project  demonstrates  graphically  how  ideas  and 
concepts  developed  in  6. 1  research  can  be  transitioned  directly  into  the  fleet  within  a  very 
short  time. 
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Abstract 

When  active  sonar  systems  are  used  to  insonate  the  sea  surface,  anomalous 
scattering  is  observ^  in  the  form  of  enhanced  backscatter,  and  more  im¬ 
portantly,  in  the  form  of  discrete,  bright  echoes.  The  most  plausible  ex¬ 
planation  for  these  effects  is  the  increased  scattering  resulting  from  the 
presence  of  bubble  plumes  and  clouds,  produced  near  the  surface  by 
breaking  waves.  This  paper  describes  some  preliminary  calculations  of 
the  backscattered  target  strengths  expected  on  the  basis  of  resonance 
scattering  from  bubble  clouds. 


I.  Introduction 

In  tests  of  low  frequency  active  sonar  systems,  false  targets  have  arisen  when  insonation 
of  the  sea  surface  is  attempted,  especially  in  circumstances  of  high  sea  state  [Gauss  et  al., 
1992].  The  origin  of  these  false  targets  is  still  unclear,  although  the  most  likely  candidates 
are  assemblages  of  gas  bubbles  in  an  acoustically  compact  form.  Bubble  clouds  are  a 
common  occurrence  in  the  near  surface  of  the  ocean  when  breaking  waves  are  present; 
these  clouds  are  likely  to  result  from  bubble  entrainment  during  wave  breaking. 
Collections  of  bubbles  in  diffuse  concentrations  in  the  form  of  plumes  have  been  ob¬ 
served  at  depths  of  several  meters  [Monahan,  1971;  Thorpe,  1982;  Farmer  and  Vagle, 
1988],  presumably  drawn  to  these  depths  as  a  result  of  convective  flows  such  as  Langmuir 
circulation  and  thermal  mixing  [Thorpe,  1982].  Significant  acoustic  backscattering  from 
the  sea  surface  can  result  either  from  these  small,  relatively  dense  clouds  that  are  near  the 
surface,  or  from  the  larger,  relatively  diffuse  plumes  that  can  extend  to  greater  depths. 

It  has  been  shown  that  the  available  surface  scattering  data  [Chapman-Harris,  1962; 
Ogden  and  Erskine,  1992]  can  be  accounted  for  in  terms  of  either  weak  scattering  (the 
well-known  Bom  Approximation)  from  deep,  diffuse  bubble  plumes  generated  by 
Langmuir  circulation  [MacDonald,  1991;  Henyey,  199J],  or  by  resonance  scattering  from 
higher  void  fraction  clouds  near  the  surface  [Prosperetti  and  Sarkar,  1992];  it  is  the  con¬ 
tention  in  this  paper  that  the  observed  large  target  strength  false  echoes  result  principally 
from  detatched  bubble  clouds;  furthermore,  we  present  in  this  paper  the  range  of  bubble 
cloud  and  environmental  parameters  that  are  likely  to  result  in  these  bright  targets. 


2.  Background 


The  attempts  to  characterize  acoustic  scattering  from  the  ocean  surface  in  the  absence  of 
bubble  clouds,  due  to  Bragg  Scattering  alone,  have  resulted  in  significant  disagreement 
between  the  calculations  [McDaniel,  1987]  and  the  experiment^  data  [Chapman  and 
Harris,  1962;  Ogden  and  Erskine,  1992].  Consequently,  MacDonald  [1991]  and  Henyey 
[1991]  have  used  weak  scattering  theory  (Bom  approximation)  to  obtain  the  surface 
backscatter  in  the  presence  of  "tenuous"  (void  fractions  less  than,  say,  10'^  %)  bubble 
plumes  of  various  configurations  and  orientations.  Their  calculations  assume  that  the 
clouds  are  sufficiently  diffuse  so  that  multiple  scattering  can  be  ignored;  consequently, 
the  scattered  sound  energy  is  mostly  specular.  Their  results  indicate  that  if  one  wishes 
only  to  account  for  the  average  surface  backscatter,  then  the  tenuous  bubble  plumes  gen¬ 
erated  by  Langmuir  circulation  are  sufficient.  However,  it  is  not  yet  clear  whether  these 
approaches  can  account  for  the  presence  of  the  "bright  echoes"  or  "hot  spots"  observed  in 
the  critical  sea  tests  and  described  by  Gauss  et  al.,  [1992].  Thorsos  [1992]  has  examined 
the  effect  of  a  rough  surface  and  not^  that  focusing  from  appropriate  contours  can  signif¬ 
icantly  enhance  the  calculated  backscatter  from  tenuous  clouds.  However,  we  shall  fol¬ 
low  a  different  approach  and  assume  that  there  are  concentrations  of  bubbles  in  the  form 
of  clouds  that  are  of  sufficient  void  fraction  to  lead  to  a  resonance  oscillation  of  the  cloud 
itself,  thus  resulting  in  high  target  strengths  at  low  frequencies. 

These  cloud  oscillations  are  called  "collective  oscillations"  and  represent  a  type  of  acous¬ 
tic  backscatter  that  is  fundamentally  different  from  that  described  by  Bom-approxima¬ 
tion,  weak-scattering  theory.  Collective  effects  occur  when  ihe  acoustic  wavelength  is 
considerably  larger  than  the  dimensions  of  the  cloud  and  the  resonance  frequency  of  the 
individual  bubbles  comprising  the  cloud  is  much  higher  than  the  insonation  frequency; 
thus,  all  the  bubbles  oscillate  essentially  in  phase.  Because  the  compressibility  of  the 
cloud  is  similar  to  that  for  a  single  gas  bubble  while  the  induced  mass  is  associated  with 
that  of  the  entire  cloud,  the  oscillation  frequency  can  be  quite  low,  and  particularly,  much 
lower  than  that  of  an  individual  gas  bubble.  An  alternative  but  equivalent  explanation  is 
that  because  the  phase  speed  in  the  bubbly  mixture  is  greatly  reduced  (sometimes  even 
below  that  for  a  pure  gas),  the  oscillation  frequency  is  correspondingly  reduced. 


2.1  Collective  Oscillations 

Carey  and  Browning,  [1988]  and  Prosperetti  [1988]  independently  suggested  that  bubble 
clouds  whose  geometrical  dimensions  were  small  with  respect  to  a  wavelength  could  be¬ 
have  as  a  compact  scatterers.  Evidence  for  the  existence  of  collective  oscillations  have 
been  firmly  established  by  laboratory  work  [Yoon,  et  al.,  1991;  Nicholas,  et  al.,  1992;  Lu, 
et  al.,  1991],  and  also  by  field  experiments  in  a  large  fresh  water  lake  [Roy,  et  al.,  1992]. 
Furthermore,  recent  data  by  Fanner  and  Ding,  [1992]  on  sources  of  ambient  noise  in  the 
ocean  provide  strong  support  for  the  existence  of  low  frequency  emissions  indicative  of 
collective  oscillations.  If  these  clouds  arc  observed  to  radiate  at  low  frequencies  by  col¬ 
lective  effects  then  it  is  likely  that  they  would  al^o  act  as  effective  scatterers  of  low  fre¬ 
quency  sound. 

2.2  Preliminary  Results 

We  have  performed  some  preliminary  experiments  of  the  low  frequency  scattering  char¬ 
acteristics  of  bubble  clouds.  The  test  plan  and  some  initial  results  are  described  in  a  pre¬ 
vious  SACLANT  Symposium  report  [Carey  and  Roy,  1993],  and  in  a  more  widely  dis¬ 
tributed  publication  [Roy,  et  al.,  1992];  they  are  shown  for  completeness  in  Fig.  1. 


Fig.  la.  Test  Plan  for  Seneca  Lake  scattering  Fig.  lb.  Low-frequency  backscatiering 

experiment.  results. 


The  principal  results  of  this  experimental  study  can  be  summarized  as  follows: 


•  Artificially  generated  bubble  clouds  of  ellipsoidal  geometry  and  about  0.5  meter 
in  diameter  and  1.0  meter  in  length  were  created  at  a  depth  of  about  90  meters  in  a  fresh 
water  lake.  The  clouds  were  insonified  with  both  a  directed  beam  (parametric  array)  and 
an  omnidirectional  conventional  source  over  a  frequency  range  from  200  Hz  to  approxi¬ 
mately  14  kHz. 

•  Measurements  of  the  target  strength  (TS)  as  a  function  of  frequency  show  rela¬ 
tive  maxima  at  approximately  0.3  and  1.3  kHz,  as  well  as  several  other  (higher)  frequen¬ 
cies.  The  amplitude  of  these  peaks  is  quite  large  and  indicative  of  resonance  effects. 

•  The  resonance  frequencies  of  the  individual  bubbles  comprising  the  clouds  are 
on  the  order  of  2-3  kHz,  and  are  so  much  larger  than  the  low  frequency  maximum  that 
these  low  frequency  peaks  are  most  likely  due  to  collective  oscillation  resonances.  If  the 
cloud  is  treat^  as  an  acoustically  compact  object  with  a  velocity  of  sound  significantly 
different  than  that  of  water,  then  the  fundamental  (monopole)  resonance  frequency  de¬ 
pends  solely  upon  its  volume,  and  its  effective  acoustic  impedance.  Using  a  modified 
.  Minnaert  formula,  a  calculated  resonance  frequency  of  about  324  Hz  can  be  obtained  for 
the  lowest  peak;  this  value  compares  favorably  with  the  measured  resonance  of  310  Hz 
[Roy,  et  al.,  1992;  Carey  and  Roy,  1993]. 

•  The  target  strength  of  these  clouds  insonified  near  resonance  is  on  the  order  of  0 
db.  Thus,  they  represent  bright  targets  and  compact  scatterers.  Calculations  of  the  target 
strength  based  on  resonance  scattering  [Roy,  et  al.,  1992;  Carey  and  Roy,  1993]  suggest 
values  on  this  order. 

The  success  of  these  preliminary  studies  has  emboldened  us  to  attempt  a  more  systematic 
and  detailed  analysis  of  low  frequency  resonance  scattering  from  near-surface  bubble 
clouds;  our  progress  along  these  lines  is  described  in  the  sections  to  follow. 


3.  Approach 


We  follow  the  approach  of  Morse  and  Ingard  [1968]  in  which  we  assume  a  plane  wave 
incident  on  a  compliant  sphere  of  radius  a^.  surrounded  by  a  continuous  medium  of  den¬ 
sity  and  sound  speed  p  and  c  respectively.  Likewise,  we  consider  the  sphere  to  be  a 
homogeneous  medium  of  density  and  sound  speed  p  and  c  respectively.  (It  should  be 
noted  that  the  subscript  'Y  "  refers  to  the  bubbly  mixture,  no\  the  pure  liquid.)  We  shall 
assume  that  the  target  is  illuminated  by  plane  waves  and  that  it  radiates  a  spherically  out¬ 
going  acoustic  wave. 

We  shall  also  assume  that  the  sphere  (bubble  cloud)  is  composed  of  many  bubbles;  thus, 
the  medium  is  dispersive  with  effective  density  and  wave  number  given  by, 

Pt=Ppnir+(^~P)p 

,  2  _  1,2  ,  Anoi^ari  .  (1) 
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where  k  =  (O  !  c  is  the  wave  number  in  the  liquid,  a  is  the  radius  of  individual  bubbles 
(considered  to  monodispersed  in  size),  T)  is  the  number  of  bubbles  per  unit  volume,  is 
the  resonance  frequency  of  the  bubbles,  /J  =  Atctio^  /  3  is  the  void  fraction,  and  b  is  the 
damping  constant  (Commander  and  Prosperetti,  1989;  Lu,  ei  al,  1990].  Foro) «  co^, 
one  can  show  that  the  real  portion  of  the  complex  phase  speed  in  the  mixture  is  given  by. 


c,  = 


yp^ 


yP. 


(2) 


Here,  Y  is  the  ratio  of  specific  heats  of  the  gas.  It  should  be  noted  that  for  (5  not  too  small 
or  P„  not  too  large,  the  low  frequency  phase  speed  in  Eq.  (2)  reduces  to  the  more  famil¬ 
iar  expression. 


Q  = 


(3) 


For  the  scattering  problem  in  an  infinite  medium  we  solve  the  Helmholtz  equation  subject 
to  the  boundary  conditions  of  continuous  pressure  and  normal  velocity  across  the  surface 
of  the  sphere. 


VV  =-/  (r).  (4) 

We  take  the  solution  to  be  a  superposition  of  incident  plane  wave  and  scattered  waves: 

p  =p  +  p  where  p  (5) 

I  /  I 

Morse  and  Ingard  used  an  integral  Green's  function  method  to  demonstrate  that  the  solu¬ 
tion  for  the  exterior  scattered  wave  in  spherical  coordinates  is  represented  by  an  expan¬ 
sion  in  Legendre  polynomials  and  spherical  Hankel  functions  with  appropriate  coeffi¬ 
cients. 


(6) 


p^{r)  =  ~A'^r{2m  +  l)(l  +  R„)P„{cosd)h„{kr) 
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— '  ^  V  ^m)^m(cOS0) 

/n=0 


where  the  asymptotic  form  is  given  in  the  far  field.  The  coefficient  satisfies  the 
boundary  conditions  and  describes  the  reflectivity  of  the  sphere  where. 


and 


Pn, 


PlCi 


is  the  specific 


admittance  of  the  surface.  In  this  study  we  make  use  of  the  limiring  form  of  given  in 
Eq.  (2)  which  is  not  a  complex  phase  speed  and  hence  does  not  take  into  account  damp¬ 
ing  from  the  bubbles  within  the  cloud. 

We  are  primarily  interested  in  the  backscattering;  hence,  we  take  6  —  k.  In  the  ffee-field, 
the  TS  is  given  by  [Urick,l967], 


TSff  =  20  log 


P, 
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Carey  and  Roy  [1993]  have  shown  that  for  small  ka^,  the  monopole  term  in  Eq.  (6)  can 
be  approximated  by 
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At  resonance,  this  leads  to 


1  'iyP^ 
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(9) 


which  is  a  modified  form  of  Minnaen's  equation  for  bubble  resonance. 

The  above  equations  allow  us  to  make  some  observations  concerning  the  interdependence 
of  some  of  the  more  important  parameters,  viz, 

•  The  effective  phase  speed  in  the  bubbly  mixture  increases  with  depth,  for  a 
fixed  volume  fraction  ]3;  likewise,  for  fixed  the  resonance  frequency  of  the  cloud, 
0)^,  also  increases  with  depth. 

•  As  P  decreases/increases,  both  and  c,  increases/decreases  (other 
parameters  held  fixed). 


fixed). 


As  decreases/increases,  co^  increases/decreases  (other  parameters  held 


Furthermore,  we  can  also  deduce  how  changes  in  some  of  these  parameters  affect  the  TS? 
At  resonance,  Eq(8)  reduces  to 
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Thus,  as  the  resonance  frequency  decreases,  or  the  void  fraction  becomes  larger,  or  the 
size  of  the  cloud  increases,  the  TS  increases. 


We  now  turn  to  a  determination  of  the  TS  for  a  variety  of  environmental  and  acoustic 
conditions.  In  our  calculations,  we  choose  to  obtain  the  resonance  frequencies  by  numer¬ 
ical  methods.  We  will  approximate  Eq.  (6)  by  truncating  the  series  at  the  m  =  0 
monopole  term,  since  for  the  frequencies  of  interest  in  this  paper,  the  wavelength  of 
sound  is  many  times  greater  than  the  bubble  cloud  radius  and  the  higher  order  terms  do 
not  contribute  to  the  lowest  resonance. 

To  solve  the  problenr.  of  scattering  from  targets  near  the  ocean  surface  we  make  use  of  the 
method  of  images;  in  Fig  2  below  a  diagram  is  shown  of  our  approach. 


Fig.  2.  Diagram  of  the  theoretical  model  used  in  the  determination  of  the  scattered  target 
strengths 

Here  we  treat  the  target  as  a  point  scatterer  a  distance  d  below  a  pressure  release  surface 
with  reflectivity  coefficient  p.  (-1<  p,  <0).  The  reflection  coefficient  describes  the  rough¬ 
ness  of  the  sea  surface  where  p.  =  - 1  corresponds  to  a  smooth  pressure  release  surface  and 
p.  =  0  corresponds  to  an  extremely  rough  surface  (effectively  equivalent  to  an  infinite 
medium);  re^istic  sea  states  fall  somewhere  in  between.  Rather  than  relate  p.  to  the  fre¬ 
quency  and  wave-height,  we  have  chosen  to  evaluate  the  TS  at  fixed  values  of  p  in  order 
to  generalize  the  analysis.  Using  the  method  of  images  one  obtains, 


(11a) 
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(11b) 

where  0.  is  the  surface  grazin"  angle,  and  p,  and  p,  are  the  magnitudes  of  the  incident 
and  scattered  fields  in  the  free-field.  The  single-primes  indicate  th^  fields  neglecting  the 
surface,  the  double-primes  denote  the  image  fields,  and  |p,L  =  l^i|  /^is  the  response  of 
the  cloud  to  the  incident  plane  wave,  where  TSV  =  lOlog/^  is  the  target  strength  com¬ 
puted  in  a  field.  After  some  algebra,  we  find  that  the  TS  for  a  bubble  cloud  near  the 
sea  surface  is  given  by 
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Clearly,  when  p.  =  0  the  expression  for  TS  yields  the  free-field  target  strength;  for  |i  = 
-1,  and  the  source  close  to  the  surface,  phase  cancellation  occurs,  and  the  source  behaves 
like  a  dipole. 


If  one  considers  the  limitipg  from  of  Eq.  (12),  for  kd«l,  one  observes  a  dipole  character¬ 
istic  to  the  scattered  field  in  which  the  scattered  pressure  scales  with  (dA)'*  sin^G,  where  X 
is  the  acoustic  wavelength.  This  equation  suggests  a  complex  acoustical  behavior  which 
can  be  briefly  summarized  as  follows.  High  void-fraction  clouds  which  generate  signifi¬ 
cant  TS’s  in  Ae  free  field  may  not  be  acoustically  important  because  these  clouds  tend  to 
reside  near  the  surface,  and  thus  are  subject  to  the  mitigating  effect  of  surface  dipole  can¬ 
cellation.  This  effect  is  exacerbated  by  the  fact  that  these  large  TS  clouds  tend  to  res¬ 
onate  at  low  frequencies  (i.e.  long  wavelengths),  and  the  proximity  of  the  cloud  to  the 
surface  is  defined  relative  tc  the  acoustic  wu  ’elength. 

If  one  considers  deeper  sc  .tterers,  then  one  is  necessarily  limited  (by  oceanographic  con¬ 
straints)  to  the  consideration  of  lower  void-fraction  clouds  that  will  not  have  as  pro¬ 
nounced  a  resonance  scattering  characteristic.  *ndeed,  clouds  in  the  deepest  portion  of 
the  bubble  layer,  (order  10  m)  are  very  tenuous  and  probably  do  not  resonate  at  all.  It 
seems  likely  that  there  are  optimum  combinations  of  cloud  depth,  cloud  characieristics 
and  frequency  that  produces  significant  backscatter  target  strengths  (larger  than,  say  -10 
db).  In  the  next  section,  uiring  the  equations  determined  above,  we  have  explored  a  va¬ 
riety  of  conditions  that  could  give  rise  to  significant  scattering  TS's. 


.4.  Results 

We  have  generated  a  ser-s  of  multidimensional'figures  to  display  the  calculated  target 
strengths  of  these  bubble  assemblages  as  a  function  of  several  relevant  parameters.  These 
figu^  are  quite  complex  and  require  studious  attention  in  order  to  fully  grasp  the  princi¬ 
pal  implications  of  the  data.  The  parameters  shown  in  these  figures  with  brief  comments 
where  appropriate  are  as  follows:  TS— the  target  strength,  defined  as  in  Eq.  12  above;  P- 
the  void  fraction,  defined  above;  d— the  depth  of  the  center  of  the  bubble  cloud  below  the 
ocean  surface;  Gg— the  grazing  angle  of  the  incident  sound  beam;  ac— bubble  cloud  radius; 

a  measure  of  me  reflectivity  of  the  surface:  fox  ^  =  -1,  the  surface  is  perfectly  reflect¬ 
ing,  for  |i  =  0,  the  surface  is  at  infinity  and  there  is  no  reflection  from  the  surface. 


Consider  Fig.  3,  which  is  a  four-dimensional  plot  of  the  target  strength  as  a  function  of 
void  fraction,  cloud  depth,  and  cloud  resonance  frequency.  For  this  case,  the  grazing  an¬ 
gle  is  shallow— 10°— and  the  cloud  radius  is  relatively  large-0.5  m.  We  anticipate  that 
this  case  would  correspond  to  the  insonation  of  a  "bubble  plume",  formed  from  the  con¬ 
vection  of  entrained  gas  bubbles  to  a  considerable  depth;  however,  this  plume  will  be 
treated  as  a  resonant,  compact  scatterer.  We  presume  that  the  large  bubbles  have  risen  by 
gravitational  forces  to  the  surface,  and  that  consequendy  the  remaining  bubbles,  and  par¬ 
ticularly  the  void  fraction,  are  both  relatively  sm^l.  The  data  shown  in  the  lower  right- 
hand-comer  of  Fig.  3  correspond  to  p  =  -1,  which  presumes  a  perfectly  reflecting  surface. 
Note  that  there  is  little  scattering  for  small  depths;  in  this  case  the  cloud  acts  as  a  dipole 
scatterer,. and  thus  for  small  grazing  angles,  the  TS  is  very  low.  However,  as  the  depth 
increases  to  a  few  meters,  even  for  void  fractions  as  low  as  10*^,  relatively  large  TS's  are 
observed.  For  example,  at  a  depth  of  5  m,  and  for  a  void  fraction  of  10"^,  the  observed 
target  strength  is  on  the  order  of  0  db  for  an  insonation  frequency  of  300  Hz. 

With  a  value  of  p  =  -1,  it  is  assumed  that  the  surface  is  perfectly  reflecting,  a  situation 
unlikely  to  be  realized  in  a  rough  sea.  The  data  shown  in  the  upper  left  hand  comer  of 
this  figure  corresponds  to  value  of  p  =  0,  or  for  a  surface  so  rough  that  no  coherent  energy 
is  reflected.  Note  that  in  this  case,  there  is  little  depth  dependence-no  phase  cancellation 
from  the  reflecting  surface— and  this  cloud  is  a  high  TS  scatterer  for  all  "depths".  Perhaps 
a  more  reasonable  representation  of  the  surface  effect  is  presented  in  the  lower  left-hand- 
comer,  calculated  for  a  value  of  p  =  -0.5,  which  thus  corresponds  to  the  intermediate 
case.  It  is  seen  from  this  figure,  that  TS's  on  the  order  of  0  db  should  be  observed  for 
plumes  with  void  fractions  above  10“^,  at  depths  between  2-6  m,  and  at  insonation  fre¬ 
quencies  on  the  order  of  300  Hz.  As  the  void  fraction  falls  below  lO'^,  the  backscattered 
intensity  rapidly  falls  in  magnitude. 

Let  us  now  compare  and  contrast  these  results  for  a  relatively  large,  low  void-fraction 
plume  with  the  case  shown  in  the  next  figure.  Shown  in  Fig.  4  are  plots  for  a  bubble 
cloud  of  radius  0.1  m,  and  a  grazing  angle  of  10°  (as  in  Fig.  3).  Consider  again  first  the 
case  for  a  perfect  reflecting  surface,  as  shown  the  lower  right-hand-comer.  'This  cloud  is 
relatively  small;  thus,  the  monopole  resonance  frequency  is  rather  high  (of  order  400  Hz 
near  the  surface).  Consequently,  it  can  produce  a  large  backscattered  TS  even  when  it  is 
within  a  few  meters  of  the  surface.  Since  we  can  presume  that  for  such  a  small  cloud,  it 
is  not  unreasonable  to  expect  a  large  void  fraction,  we  shall  consider  values  of  P  as  high 
as  10'3.  Note  that  for  the  surface  conditions  of  p  =  -1,  TS's  of  0  db  can  be  expected  only 
for  clouds  3-4  meters  below  the  surface. 

Consider  next,  however,  the  conditions  demonstrated  in  the  upper  left-hand-comer,  in 
which  the  surface  condition  corresponds  to  a  very  rough  surface— a  more  likely  occur¬ 
rence  for  large  void  fraction  clouds.  For  this  case,  TS's  of  0  db  can  be  observed  at  any 
depth  for  a  range  of  frequencies  near  400  Hz.  Again,  if  the  more  realistic  case  of  p  =  -0.5 
is  considered  (lower  left),  a  TS  of  0  db  can  be  expected  to  occur  for  a  cloud  with  a  void 
fraction  of  10'^,  a  radius  of  0. 1  m,  a  frequency  of  400  Hz,  and  at  a  depth  of  about  4  m. 

Finally,  let  us  consider  Fig.  5,  which  summarizes  the  principal  features  of  this  repon. 
Our  field  measurements  at  Lake  Seneca  (Roy,  et  aJ.,  1992;  Carey  and  Roy,  1993]  demon¬ 
strated  that  TS's  of  approximately  0  db  could  be  achieved  for  resonance  oscillations  of  a 
bubble  cloud,  far  removed  from  the  surface,  with  a  void  fraction  of  about  10"^,  a  radius  of 
about  0.25  m  and  at  an  insonation  frequency  of  about  300  Hz.  Funhermore,  Prosperetti 
and  his  students  [Prosperetti,  et  al.,  1993;  Sarkar  and  Prosperetti,  1993]  have  demonstrated 
that  bubble  clouds  of  similar  size  and  void  fraction,  when  located  very  near  the  surface, 
could  account  for  the  Chapman-Harris,  Ogden-Erskine  surface  acoustic  backscatter  when 


perfectly  reflecting;  when  p.  is  0,  no  reflection  occurs. 


Fig.  4.  Three  dimensional  plot  of  the  calculated  maximum  target  strength  as  a  function  of  depth  and  volume  fraction  for  a  grazing  angle  of  10^  and  a 
cloud  radius  of  0.1  m.  The  frequency  of  the  monopole  resonance  is  shown  via  the  color  bar  to  the  right.  The  effect  of  the  surface  is  described  via  a 
reflection  coefTicient,  tu  When  \i  is  -1,  the  surface  is  perfectly  reflecting;  when  \x  is  0,  no  reflection  occurs. 


Backscatter  TS  at  Monopole  Resonance 
Radius  =  O.lm  Beta  =  0.005 

(350-500)  Hz 


treated  as  resonant  scatterers.  We  show  now  in  Hg.  5  the  expected  individual  resonance 
scattering  characteristics  of  these  clouds  as  a  function  of  such  parameters  as  grazing 
angle,  depth  below  the  surface,  and  surface  conditions. 

Shown  in  the  lower  right-hand  comer  of  Fig.  5  is  the  predicted  backscattered  target 
strength  for  a  bubble  cloud  at  its  monopole  resonance  (varying  between  350-500  Hz),  for 
a  radius  of  0.1  m,  a  void  fraction  of  0.5  %,  and  at  a  perfectly  reflecting  surface.  Also 
shown  in  this  figure  are  the  cases  for  a  free  field  (upper  left)  and  for  the  intermediate  case 
of  a  partially  reflecting  surface  (lower  left).  Consider  the  plot  in  the  lower  left.  This  case 
represents  our  estimate  of  the  conditions  expected  to  give  rise  to  false  echoes  in  low-fire- 
quency,  near-surface  scattering.  Note  that  one  should  expect  echoes  with  TS's  of  0  db  for 
a  wide  range  of  grazing  angles  and  cloud  depths.  Bright  echoes  (of  order  0  db)  can  be 
expected  to  occur  for  grazing  angles  from  10-70°,  and  for  cloud  depths  from  near  the  sur¬ 
face  to  10  meters.  This  figure  indicates  that  the  ocean  surface  is  very  rich  with  possibili¬ 
ties  for  bright  echoes  from  resonant  bubble  clouds. 

The  conditions  that  lead  to  the  production  of  bubble  clouds  (high  winds  and  high  sea 
states)  also  favor  the  generation  of  high  scattering  TS's.  First  of  all,  high  sea  state  means 
that  the  ocean  surface  will  be  rough;  thus,  |i=-0.5  is  not  an  unrealistic  approximation;  sur¬ 
face  cancellation  won't  be  as  important  except  for  very  low  frequencies.  Secondly, 
stormy  conditions  usually  imply  an  unstable  water  column,  which  is  often  upward  refract¬ 
ing.  TTius,  the  incident  angles  become  steeper  and  the  sin^G  dependence  is  no  longer  as 
important  These  two  effects  tend  to  make  far-field-like  behavior  much  more  pronounced. 

Finally,  we  make  note  of  the  experiments  of  Lamarre  and  Melville  [1992]  who  obtained 
measurements  of  the  void  fractions  of  ocean-generated  bubble  clouds  in  the  open  ocean 
off  the  coast  of  Delaware  during  the  last  week  of  February,  1991.  Their  data  indicate  that 
significant  numbers  of  clouds  are  produced  with  void  fractions  on  the  order  of  0.5  %; 
they  remark  that  their  data  are  consistent  with  their  laboratory  results  and  are  several  or¬ 
ders  of  magnitude  higher  than  the  often-reported  time-averaged  values  [Farmer  and 
Vagle,  1988]. 

Further  experiments  are  required,  however,  to  determine  the  precise  acoustic  and  envi¬ 
ronmental  conditions  that  give  rise  to  these  false  echoes. 

5.0  Summary 

False  targets  can  occur  for  low  frequency  sound  scattering  from  the  sea  surface.  It  is  un¬ 
likely  that  Bragg  scattering  from  a  rough  surface  is  the  source  of  these  echoes;  rather, 
bubble  clouds  or  plumes  resulting  from  breaking  waves  represent  a  more  plausible  expla¬ 
nation.  We  have  shown  that  bubble  clouds  or  plumes  can  act  as  resonant  scatterers  at  low 
frequencies  and  generate  target  strengths  as  large  as  0  db  for  a  wide  variety  of  conditions 
that  arc  expected  to  occur  in  rough  seas  with  wind-driven  breaking  waves.  Future  exper¬ 
iments  are  planned  to  determine  more  precisely  th^  measured  targets  strengths  in  terms  of 
the  environmental  conditions  that  would  generate  bubbles  clouds  of  the  required  size, 
shape  and  void  fraction  to  produce  these  false  targets. 
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Acoustically  Active  Surfaces 

Acoustically  Active  Sutfaces  and  Active  Noise  and  Vibration  Control 
Principal  Investigator:  F.  Douglas  Shields 


Research  Accomplished  in  1992: 

This  research  began  in  1988  with  an  effort  to  develop  smart  acoustically  active  coatings  using 
piezorubber  and/or  PVDF.  Such  smart  surface  coatings  are  constructed  by  combining  a  sensing 
layer  with  one  or  more  actuating  layers.  The  digital  signal  processing  capability  has  been 
developed  that  enables  the  signal  from  a  sensing  transducing  layer  to  be  passed  through  a  digital 
signal  processor  and  applied  to  one  or  two  actuating  layers.  With  one  actuator,  it  has  been 
possible  to  control  either  transmission  or  reflection  of  plane  waves.  With  two  actuating  layers  it  is 
possible  to  simultaneously  control  both  reflection  and  transmission. 

The  work  with  smart  coatings  has  expanded  to  include  three  aspects.  First  concerns  the 
development  of  actuators  and  sensors,  especially  those  for  use  at  low  frequencies;  the  second 
concerns  the  study  of  different  control  algorithms;  and  the  third  involves  the  study  of  the  coupling 
between  the  active  surface  and  the  medium.  This  coupling  can  be  complicated  when  the  velocity 
of  the  surface  becomes  position  dependent  The  work  is  these  three  areas  is  discussed  below. 

The  digital  signal  processors  needed  for  the  "sman"  surface  control  are  the  same  as  those 
needed  for  active  noise  reducing  headsets.  About  three  years  ago  a  research  project  was 
undertaken  with  the  Navy  Experimental  Diving  Unit  in  Panama  City,  Rorida  to  develop  such  a 
head  set  for  divers  to  use  to  reduce  noise  in  a  diving  helmet. 

A  prototype  headset  with  its  digital  control  was  delivered  to  the  NEDU  in  February.  Work 
was  begun  in  February  on  an  SBIR  with  the  Army  to  develop  an  ANR  stethoscope  for  use  in 
helicopters  and  other  noisy  vehicles.  This  is  also  now  in  the  prototype  stage. 


n.  Summary  of  Work 

Progress  in  the  three  areas  of  research  with  smart  coatings  funded  by  ONR  is  discussed  in 
this  section. 

Sensor  Development 

Because  of  their  limited  thickness,  active  coatings  are  intrinsically  insensitive  at  low 
frequencies.  In  1991  NCPA  developed  panels  consisting  of  transflexural  capsules  potted  in  a 
matrix  of  polyurethane.  During  the  interval  covered  by  this  report,  these  panels  were  used  to 
eliminate  reflections  from  a  pressure  release  surface  in  a  water  filled  pulse  echo  tube  10  inches  in 
diameter  and  20  feet  long  at  frequencies  as  low  as  100  Hz.  In  order  to  predict  the  performance  of 
such  panels,  the  three  port  impedance  matrix  of  these  panels  was  measured  experimentally  using 
a  Laser  Doppler  Velocimeter.  Work  was  begun  on  a  contract  with  The  Sandia  National 
Laboratory  to  collaborate  in  the  development  of  transflexural  capsules  that  could  be  used  at  great 
depths  in  the  ocean.  As  part  of  that  contract,  the  performance  of  a  number  of  different  capsules 
was  measured  as  a  function  of  pressure.  An  analytic  expression  was  developed  to  calculate  their 
performance  and  preliminary  work  was  done  on  a  finite  element  analysis  of  the  unit. 

Algorithm  Development 

The  control  of  smart  coatings  at  pressure  release  surfaces  is  especially  susceptible  to 
instability  due  to  degenerative  feedback  from'a  pressure  sensor.  With  a  velocity  sensor  the 
feedback  is  nondegenerative.  The  pressure  and  velocity  sensor  reverse  rolls  at  a  rigid  surface. 
This  is  one  reason  proposed  for  using  a  velocity  sensor  with  a  soft  surface.  However,  a  velocity 
sensor  that  does  not  alter  the  surface  and  has  the  needed  sensitivity  is  very  difficult  to  configure. 
During  the  period  covered  by  this  report,  an  analysis  of  the  instability  problems  associated  with 
various  control  algorithms  was  made  and  the  results  displayed  on  three  dimensional  plots  of  the 
reflection  and  transmission  coefficients  as  a  function  of  complex  gain.  A  paper  was  prepared 
discussing  these  caleulations  and  the  results  of  some  experimental  work  with  multilayered 


surfaces.  This  paper  was  ruled  by  ONR  to  contain  sensitive  material  and  it  was,  therefore,  not 
published. 

In  the  area  of  active  noise  control,  innovative  work  was  done  in  developing  both 
frequency  domain  and  time  domain  adaptive  filters.  Some  of  this  work  was  described  in  a  paper 
that  was  submitted  for  publication  to  the  Journal  of  the  Acoustical  Society  of  America,  but  is  still 
in  the  review  process. 

Coupling  Between  Active  Surfaces  and  the  Medium 

To  control  reflections  of  plane  waves  that  strike  a  plane  surface  obliquely  requires  a 
phased  array  of  active  surface  elements.  Tim  Ruppel  built  such  a  surface  as  a  part  of  his 
dissertation  project.  He  made  a  plane  surface  with  active  strips  and  was  able  to  drive  these  strips 
with  a  signal  whose  phase  increased  at  a  controlled  rate.  Using  this  surface  he  controlled 
reflections  of  plane  waves  ft'om  the  surface  and  radiated  plane  waves  at  specified  angles.  This 
work  was  published  and  is  referenced  below.  Additional  work  has  been  initiated  to  construct  a 
square  wave  guide  with  inner  walls  that  can  be  driven  with  a  controllable  phase  shift  so  as  to 
propagate  waves  with  a  controllable  frequency  and  wave  number  on  iilflt'i  WUAl. 

A  study  has  also  been  initiated  by  Dr.  Lafleur  of  the  propagation  of  different  modes  in  a 
fluid  filled  tube  with  elastic  walls. 
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ABSTRACT 


A  stochastic  scattering  model  is  derived  in  which  the  backscatter  from  the  near- 
surface  oceanic  bubble  layer  is  written  directly  in  terms  of  an  integral  over  the  wavenumber 
spectrum  ("power"  spectrum)  of  the  sound  speed  fluctuations  in  the  layer.  A  factored  form 
is  given  for  the  integral  that  allows  the  backscatter  cross  section  per  unit  area  to  be 
expressed  as  the  product  of  a  "geometric  factor"  and  an  effective  horizontal  wavenumber 
spectrum.  Because  a  power  spectrum  formulation  is  statistical,  there  are  no  assumptions 
about  the  geometry  of  the  bubble  layer.  (E.g.,  hemispherical  or  cylindrical  plumes  are  not 
assumed.)  By  dividing  measured  data  for  backscatter-vs-frequency  by  the  geometric 
factor,  we  inverted  scattering  data  from  12  different  deep-ocean  reverberation 
measurements  and  directly  inferred  the  wavenumber  spectrum  of  the  sound  speed  in  the 
bubble  layer.  For  all  12  measurements,  the  inferred  wavenumber  spectrum  is  an  inverse 
power  law  of  the  form  P(K)=AK'P  where  A  is  a  strength  parameter,  K  is  the  horizontal 
wavenumber  for  a  Fourier  component  of  the  sound-speed  distribution,  and  the  mean  value 
of  the  spectral  roll-off  exponent  p  is  3.86  ±  .45  .  The  consistency  in  the  inferred 
wavenumber  spectrum  strongly  suggests  that  on  scales  of  less  than  half  an  acoustic 
wavelength  (5  to  10  m),  the  sound-speed  structure  in  the  bubble  layer  is  governed  by 
turbulence  in  the  inertial  subrange  (Kolmogorov  subrange)  which  has  a  universal  value  of 
P  =  1 1/3  =  3.67  for  fully  developed  isotropic  turbulence.  Using  Von  Karman’s 
inteipolation  formula  for  Kolmogorov  turbulence,  together  with  occanogrju)hically 
constrained  input  parameters,  we  computed  theoretical  backscatter  cross  sections  and 
compared  them  with  the  empirical  fit  of  Ogden  and  Erskine  to  Critical  Sea  Test  data.  It  is 
shown  that  with  no  adjustable  parameters,  the  stochastic  scattering  model  gives  a  good 
account  of  the  observed  backscatter,  as  a  function  of  both  frequency  and  grazing  angle. 


INTRODUCTION 

In  the  scattering  of  sound  from  the  sea  surface,  a  significant  discrepancy  exists 
between  the  observed  measurements  of  backscatter  and  the  values  predicted  by  the  classic 
theory  of  Bragg  scattering  from  rough  sea  surfaces .  At  high  sea  states  and  low  grazing 
angles,  observed  backscattering  strengths  are  10  to  100  times  stronger  than  the  strengths 
predicted  by  Bragg  scattering  theory  applied  to  sea  surfaces  having  realistic  scales  of 
roughness  (McDaniel,  1988, 1993;  Ogden  and  Erskine,  1992).  It  is  doubtful  that  the 
discrepancy  can  be  removed  by  improvements  to  rough  surface  scattering  theory  or  by 
adjustments  to  the  known  sea  surface  roughness  spectra.  Consequently  it  has  been 
hypothesized  that  the  observed  high  reverberation  levels  are  due  not  to  rough  surface 
scattering  but  to  scattering  from  the  near-surface  oceanic  bubble  layer  which  is  created 
when  air  is  entrained  by  breaking  waves  (Urick  and  Hoover,1956;  McDaniel,  1988,  1993). 

Since  large  bubbles  rapidly  rise  to  the  surface,  the  bubble  layer  itself  is  composed 
primarily  of  clouds  of  very  small  bubbles  ("microbubbles")  that  are  distributed  and  mixed 
by  downwelling  currents  and  turbulence  (Thorpe,  1982, 1984;  Thorpe  and  Hall,  1983; 
Osborn,  et  al.,1992).  The  microbubbles  are  transported  essentially  passively  until  they 
finally  go  into  solution,  usually  after  many  minutes  (Thorpe,  1982;  Thorpe  and  Hall, 

1983).  Although  the  microbubble  layer  has  a  low  density  of  bubbles  (air  volume  fractions 
of  10'^  or  smaller),  it  nevertheless  has  a  marked  effect  acoustically  (Farmer  and  Vagle, 
1989;  Su,1992;  CartmiU  and  Su,  1992;  Su  et  al.,  1993). 

At  low  frequencies  where  the  scattering  is  nonresonant  (f  <  5  kHz),  the  sound 

speed  in  the  bubble  layer  can  be  considered  to  be  the  that  of  a  water-bubble  mixture.  Small 
concentrations  of  bubbles  significandy  increase  the  compressibility  of  a  water  parcel 
without  measurably  changing  its  density.  Since  sound  speed  is  inversely  proportional  to  the 
square  root  of  density  times  compressibilty,  the  sound  speed  in  the  bubble  layer  is 
decreased.  For  example,  for  a  volume  fraction  of  lO'^,  the  speed  of  sound  is  reduced  by 
about  170  m/s. 
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On  records  from  upward-looking  sonars,  the  distribution  of  bubbles  in  the  layer 
appears  to  be  random  in  space  and  time.  (Thoipe,1982;  Farmer  and  Vagle.  1989;  Farmer, 

1992;  Commander,  1992).  Thus  the  acoustic  "picture"  of  the  near-surface  bubble  layer  is  a 
region  with  a  stochastic  distribution  of  sound  speed  that  is  determined  by  upper  ocean 
dynamic  processes. 

In  this  article  we  propose  that  the  clouds  of  microbubbles  discussed  above  are  an 
important,  and  perhaps  the  primary,  physical  mechanism  that  controls  the  backscatter  at 
low  grazing  angles  and  high  sea  states.  Such  a  proposal  has  been  made  already  by  other 
workers  (McDonald,  1991;  Henyey  1991).  In  fact,  the  work  described  here  was  directly 
motivated  by  the  original  research  of  McDonald.  Unlike  previous  workers,  however,  we 
do  not  assume  the  bubble  layer  to  be  a  collection  of  plumes  with  a  specified  geometry  such 
as  cylinders  or  hemispheres.  Rather,  we  more  realistically  consider  the  bubble  layer  to  be  a 
stochastic  medium  with  certain  statistical  properties.  With  this  picture  of  the  bubble  layer, 
the  backscatter  is  viewed  as  Bragg  scattering  from  the  Fourier  components  of  the  sound- 
speed  fluctuations  created  by  the  stochastic  microbubble  distribution.  As  is  discussed 
below,  this  approach  makes  it  possible  to  relate  the  backscatter  directly  to  the  wavenumber 
spectrum  of  the  sound-speed  fluctuations  in  the  bubble  layer,  thereby  connecting  the 
backscatter  to  a  realistic  oceanogrt^hic  process  -  turbulence. 


1.  TOEOREnCAL  FORMULATION 

We  want  to  calculate  the  backscatter  from  the  near-surface  microbubble  layer  under 
the  assumption  that  the  bubble  distribution  varies  stochastically  in  space  and  time.  To 
arrive  at  such  a  formulation  we  represent  the  bubble  layer  as  an  effective  fluid  continuum 
with  a  fluctuating  sound  speed  that  on  the  average  is  less  than  the  sound  speed  in  bubble- 
free  water.  We  first  outline  the  basic  assumptions  of  the  bubble-layer  model.  Then  an 
expression  is  derived  for  the  backscatter  one  obtains  fix)m  the  bubble  layer  model.  Finally, 
we  approximate  the  backscatter  cross-section  per  unit  area  as  the  product  of  a  "geometric 
factor"  times  the  effective  horizontal  wavenumber  spectrum  for  sound  speed  fluctuations  in 
the  bubble  layer, 

A.  A  Stochastic  Bubble-Layer  Model 

When  a  wave  breaks  in  the  ocean,  a  wide  spectrum  of  bubble  sizes  is  generated. 
The  larger  bubbles  rise  to  the  surface  rt^idly,  leaving  behind  a  turbulent,  smoke-like 
cloud  of  microbubbles  with  radii  of  10-100  microns  (Su  and  Cartmill,  1992;  Su  et  al., 
1993).  Since  the  microbubbles  may  stay  in  the  water  for  some  time  before  dissolving,  they 
can  be  transported  downward  by  turbulent  downwelling  currents  (Thorpe,  1982,  1984; 
Thorpe  and  Hall,  1983;  Osborn,  et  al.,  1992).  Because  of  turbulence,  the  microbubble 
layer  is  a  stochastic,  inhomogeneous,  time-dependent  layer  and,  in  high  sea  states,  can 
extend  as  far  as  20  m  below  the  sea  surface  (Farmer,  1992). 

As  discussed  in  the  Introduction,  the  effect  of  the  microbubble  layer  on  low- 
frequency  sound  is  to  "soften"  the  water  so  that  the  sound  speed  is  lowered  and  the 
wavenumber  is  increased.  Thus,  in  a  continuum  nrodel,  the  square  of  the  total 
wavenumber,  Kli  ,  that  appears  in  the  wave  equation  can  be  written  as  the  contribution 
from  bubble-free  water  plus  the  contribution  from  the  bubble  distribution: 


(1) 


6 


vv^f^re  ^  oj^/c^  ic  tho  wavCiiU.nbCi  fcr  bubbL-free  seawater,  anil  kg  is  uie  change  due  to 
bubbles.  In  terms  of  the  bubble  number  density,  n,  the  contribution  of  the  microbubbles  to 
the  total  wavenumber  is 


k|  =  4k<  t  >  n  (2) 

where  <  t  >  is  the  average  bubble  scattering  strength  (Foldy,  1945;  Carstensen  and  Foldy, 
1947;  Morse  and  Feshbach,  1953).  At  the  frequencies  we  consider  the  micobubbles  are 
nonresonant,  and  Eq.(2)  reduces  to  the  wavenumber  one  obtains  using  the  sound-speed 
from  simple  mixture  theory.  Hence,  for  a  small  sound-speed  change  Ac  due  to  bubbles, 
k|  can  be  written  as 

ki  =  -2k2(Ac/co)  (3) 

where  simple  mixture  theory  gives  Ac/co  =  -  11,500^  ,  for  a  void  fraction  ^  much  less 
than  one,  (Urick,1975).  In  this  article  we  keep  in  mind  that  k|  is  due  to  a  distribution  of 

discrete  bubbles  but  always  use  a  continuum  model  to  represent  the  bubble  layer  and  to 
calculate  backscatter. 

Because  the  microbubbles  are  distributed  by  turbulent  diffusion,  at  any  given  time 
the  bubble  distribution  is  quite  inhomogeneous.  Although  experimental  measurements 
indicate  that  the  average  concentration  of  bubbles  decays  approximately  exponentially  with 
depth  ,  it  is  clear  from  the  uplooking  sonar  records  (Thorpe,  1982;  Farmer  and  Vagle, 
1989;  Farmer,  1992)  and  in-situ  measurements  (Su,  1992;  Cartmill  and  Su,  1992,  Su  et 
al.,  1993)  that  there  are  significant  fluctuations  about  the  average,  both  horizontally  and 
vertically.  Hence  we  write  k|  as  the  sum  of  an  average  value  and  a  fluctuating  part  that 
dqjends  stochastically  on  both  depth  and  horizontal  distance. 

ki  =  <ki(2)>  +6k|(r) 


(4) 
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where  <k|(z)>  =  -  2  [<Ac(z)>/co)]  is  the  average  profile  and  Sk|(r)  is  the  fluctuation 

about  Because  cf  b”bb!e  dissolution,  both  <k|(z)>  and  6k^(r)  axe  cbbeivcd  to 

decrease  with  increasing  depth  (Su,1992;  Cartmill  and  Su,  1992;  Su  et  al.,  1993). 

For  sound-speed  fluctuations  that  are  small  relative  to  Co  we  can  write  5k|(r)  as 

Sk2(r)  =  -  2k2  [Sc(r)/co]  (5) 

where  Sc(r)  is  the  fluctuation  about  the  average  sound  speed.  To  account  for  the  decrease  in 
the  bubble  concentration  with  depth  and  the  coiresponding  decrease  in  the  size  of 
fluctuations,  we  write  the  sound-speed  fluctuation  as  the  product  of  a  monotonically 
decaying  reference  function  ff»f(z)  and  a  fluctuation  parameter  t(r):  8c(r)  =  frof<z)c(r).  The 
reference  function  is  chosen  so  that  the  standard  deviation  in  sound  speed  fits  measured 
values.  We  therefore  let  fnf(z)  =  <Jc(z)i  where  Oc(z)  is  the  measured  rms  sound-speed 
fluctuation  with  depth.  The  fluctuation  parameter  t(r)  is  taken  to  represent  an  isotropic 
stationary  stochastic  process  where  the  autocorrelation  function  <c(r)t(r’)>  =  C(|^^)  is 
normalized  so  that  C(0)  =  1 .  Conceptually,  c(r)  accounts  for  the  statistics  of  the  fluid 
motion  that  transports  the  bubbles.  As  discussed  later,  the  fluid  motion  statistics  are 
directly  related  to  the  distribution  statistics  for  "conservative  passive  additives"  (Taitarski, 
1961)  that  move  with  the  fluid  as  infinite-lifetime  tracer  particles.  The  scaling  function  <Tc(z) 
accounts  for  both  the  surface  bubble  concentration  and  the  decay  with  depth  due  to  the 
finite  lifetimes  of  the  microbubbles.  Also  any  buoyancy  effects  are  included  implicitly  in 
an  empirically  determined  scaling  function. 

With  the  above  formulation,  we  have 


Skg(r)  =  -  2k^  [oc(z)/co]c(i:) 


(6) 


In  the  calculations  presented  later  we  wUl  fit  the  measured  values  of  adz)  with  an 
exponential  so  tliat  adz)  -  -  oc(0)o^p(-2/L),  where  Oc(0)  is  the  rms  sound  spcca 
fluctuation  at  the  surface  and  L  is  the  e-folding  distance.  Thus,  with  the  above  exponential 
scaling  of  the  fluctuations,  we  have  finally 

kwL,=  +  <ki(2)>  -  2ki  [  ac(0)/co  ]  e-*/L  e(r)  (7) 

The  above  expression  for  the  total  wavenumber  is  used  in  a  wave  equation  to  describe 
backscattering  from  the  oceanic  bubble  layer. 

B.  Solution  for  the  Backscattered  Field 

The  wave  equation  to  be  solved  for  the  acoustic  pressure  is  given  by 

+  Kii  0  <*) 

where  the  total  field  T  is  the  sum  of  the  field  in  the  absence  of  a  bubble  layer,  'F©  ,  and  the 
field,  ,  scattered  by  the  layer,  i.e.,  'i'  =  To  + 

The  integral  equation  associated  with  the  differential  wave  equation  is  the 
Lippmann-Schwinger  equation  (Rodberg  and  Thaler,  1967). 

i'(r)  =  i'oW  +  ^ Joofrf)  kVt)  i'(r')<f>r  (9) 

where  the  Green's  function  Oq  is  the  point  source  solution  in  the  absence  of  a  bubble  layer 
and  kKr*)  denotes  <k|(z')>  +  Sl^(r').  The  average  profile  <k|(z')>  does  not  contribute  to 
backscatter.  Consequently,  the  backscattered  field  T ,  is  given  by 


(10) 
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hor  isospec^  bubble-finee  water  and  a  pressure-release  boundary  condition  at  a  smooth  sea 
surface  we  have  for  Gq: 


Oo(i;.r;)  = 


R+ 


The  quantities  R+  and  R-,  are  defined  as 


(11) 


R±=ir-r±i  (12) 

where  with  unit  Cartesian  vectors  (1  IS  )  we  have  r±  =  fh  ±  z  and  =  x  1  +  y  ^ . 

^  ^ 

Note  that  the  vector  th  is  the  horizontal  componait  of  the  radius  vector. 

We  are  interested  in  the  far  field  so  we  approximate  Oq  as 

Oo(r,r')  =  -  ^  2i  #m  (kV  •  ''h  (13) 


where  kV  =  kgsinO  and  k'h  =  kgCosG  iS'h  are,  respectively,  the  vertical  wavenumber 
and  vector  horizontal  wavenumber  of  the  backscattered  field. 

For  weak  scattering  from  the  bubble  layer,  the  backscattered  field  makes  only  a 
small  contribution  to  the  total  field.  Hence  in  the  integral  in  Eq.(lO)  we  can  use  the  first- 
order  Bom  approximation,  which  approximates  the  total  field  in  the  bubble  layer  as  the 
unscattered  field,  i.c.,  T'  »  fo-  Taking  i'o  to  be  a  plane  wave  perfectly  reflected  from  the 
surface  we  have 

'i'o  = 


=  2i  sin  (ky  z)  •  'h 


(14) 


where  ky  =  -kgsinG  and  k^  =  kgCosG  l^h  respectively,  the  vertical  wavenumber  and 
vector  horizontal  wavenumber  of  the  unscattered  field.  With  the  results  from  Eqs.(13)  and 
(14),  the  Bom  approximation  for  the  backscattered  field  ^.(r)  is  given  by 

^.(r)=  ^  Jsm(k*vZ-)sin(kvz')^k2B(OeWh  ';h  dV  (15) 

where  k'y  =  -  ky  ,  and  qh  =  kh  -  k  h  =  2kh  . 

^  ^  ^ 

The  ensemble  averaged  backscattering  cross-section  is  given  by, 

<t  =  <|f«:atP>  (16) 

where  the  scattering  amplitude  fjcju  is  the  coefficient  of  the  spherical  wave  expfik^  r )/  r  in 
Eq.(15): 

^  Jsin2(kyz')  8k2B(r‘)  ei^h '  'h  dV  (17) 


Hence  for  backscattering  we  have, 
_  .  -tt^oKo) .  f  f 


=  [  V  '  ]  f  J  <' h  - '■h>g(kv.O**(k,,z  )C(r'-r’)  <f>r"  <f>r'  (18) 


where  g(ky,z)  s  tin2(kyz)e~*/^  and,  as  discussed  earlier,  C(r"-r')  =  <E(r'')c(r‘)>  is  the 

0*0*  00  0j 

autocorrelation  function  for  the  stochastic  quantity  t(r)  which  represents  the  effects  of 

0* 

turbulence  and  controls  the  variability  in  kg.  We  now  define  r"-  f'  as  S  =  (Sh,  Sy)  and 

00  00  0*  0S0 

rewrite  Eq.  (18)  as 


o  =  [- 


!?1(0) 


]  f  I  •‘’h-  <'"li  -  V*(kv.z'+  S,)g*(kv,z')C(Sh,Sv)rf>S  .fth'dz'  (19) 

*  *  fc/  0*  ev  ^ 


The  integral  over  z'  can  be  written  as 
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J  g(kv,z'+  Sv)g*(kv.2*)dz-  =  ^  J  |g(kv.qv)P  €‘%^dqY  (20) 

where  g(kv,qY)  is  the  Fourier  transform  of  g(kY,  z).  Thus  a  becomes 

^  ^  ]  J  |g(kY.qY)P  6(q)  dqY  dV  (21) 

where  the  normalized  wavenumber  spectrum  <5(q)  is  the  Fourier  transform  of  the 
normalized  autocorrelation  function  ,  C(S) : 

<i(q)=  fc(S)eiq  Sd3S  .  (22) 

and  q  =  (qh*  qv).  The  full  wavenumber  spectrum  for  sound-speed  fluctuations  in  the 

yv  ^ 

bubble  layer  is  proportional  to  the  normalized  spectrum. 

Experimentally,  the  quantity  measured  is  the  the  cross  section  per  unit  area,  S  , 
where  is  related  to  o  by 

o  =  J  a  d2ih'  (23) 

From  Eqs,  (23)  and  (21)  we  can  see  that  6  is  given  by 

^  ^  5  J  lt(J^v,qv)P  <i(q)  dqY  (24) 

We  want  to  use  Eq.(24)  and  surface  backscatter  measurements  to  estimate  the  horizontal 
wavenumber  spectrum  (i.e.,  power  spectrum)  of  sound-speed  fluctuations  in  the  bubble 
layer.  Consequently,  we  would  like  to  have  an  expression  for  S  that  is  factored  into  the 
product  of  a  "geometric  factor”  times  the  power  spectrum.  To  obtain  a  factored  form  we 
must  make  some  approximations  to  the  integral  in  Eq.(24). 

C.  Factored  Form  for  the  Scattering  Cross  Section 


As  indicated  in  Eq.(24),  computing  <$  requires  a  weighted  integration  over  <^(q) 
with  lfi'(kv,qv)p  as  the  weighting  function.  We  assume  that  the  fluctuations  are  separately 
isotropic  horizontally  and  isotropic  vertically  so  that  we  can  write  <2(q)  as  a  function  of  q2 
andq^:  C(q)  =  <^(q2i  q^)- of  integration,  C(q)  is  a  smooth  function 
of  q2  (such  as  a  power  law),  then  in  Eq.(24)  we  can  expand  ^(q)  in  a  Taylor  series  in  q2  . 
We  expand  0(q^,  q^)  about  some  as  yet  unspecified  average  value  <q^  : 


!)  =  CdS,  <q?>)  + 1|  ( 15  -  <i5>)  +  i 


( q^  -  <q3>)2  + 


(25) 


The  second  term  vanishes  if  we  define  <q2>  to  be  the  average  value  of  q^  over  the 
integration  interval :  <q2>  =  I1/I2,  where 


00 

=  jqr  lg(kv,qv)pdqv  =  -2n  fg(kv,z)  ^  g(kv,2)d2 
—00  J 


—00 

00 


*2  =  J  lg(kv,qY)p4qv  =  2«  J  |g(kv,2)pd2 


(26) 


(27) 


Evaluating  the  integrals  Ii  and  I2  and  simplifying  yields. 


<q?>  =  ^  ( 1  +  4 1^  l2) 


(28) 


Keeping  the  first  two  terms  in  the  Taylor  series  expansion  and  evaluating  the  integral  over 
qv,  gives  us  an  approximation  for  <9  that  is  in  factored  form: 


1 

2n 


][  J  lj(kv,<iv)P  Aiv  ]  <=(,2,  <,2>) 

n  Co 


(29) 


<^(<1^  <q3>) 


^  3L  Kam 

4ft2  L  c|  -'(1+  +  k2L2) 

where  the  result  in  Eq.(27)  has  been  used  to  evaluate  the  integral  over  qy  Note  that  the 
approximation  in  Eq.(29)  is  equivalent  to  zqDproximating  the  average  value  of  a  function 
f(x)  as  <f>  =  f  (<x>).  We  have  tested  the  approximation  in  Eq.(29)  by  comparing  with  an 
exact  numerical  evaluation  of  the  integral.  For  small  grazing  angles  and  a  power-law 
spectrum  of  the  form  <^^(qh,q?)  =  A  /(Bqg  +  Cq2)P/2,  where  A,  B,  and  C  are  constants, 
the  approximation  agrees  very  well  with  an  exact  numerical  evaluation  of  the  integral. 

We  now  define  the  effective  horizontal  wavenumber  spectrum  for  sound-speed 
fluctuations  in  the  bubble  layer  as 

P(qh)  =  <yl{0)Ci<il  <q2>)  (30) 

where  o2(0)  is  the  rms  sound-speed  fluctuation  at  z=0.  We  further  define  a  "geometric 
factor"  F 


p,  ^ 1 

c2  (1+  41^L2)(1  +  1i2l2)  ^ 


(31) 


With  the  P  and  F  defined  as  given  above,  we  have 


S  =  F  P(qh)  (32) 

Very  little  is  presently  known  about  the  wavenumber  spectrum  for  sound-speed 
fluctuations  in  the  oceanic  bubble  layer.  In  the  next  section  we  use  Eq.(32)  to  "invert" 
measured  backscatter  data  and  obtain  empirical  results  for  P(qh)- 


n.  THE  WAVENUMBEr  SPECTRUM  FOR  SOUND-SPEED  FLUCTUATIONS  IN 
THE  BUBBLE  LAVEx^ 

The  stochastic  model  discussed  above  has  been  used  to  invert  surface  backscatter 
data  from  12  different  deep-ocean  measurements.  As  will  be  shown,  all  the  measurements 
yield  an  inverse  power-law  wavenumber  spectrum  that  is  consistent  with  turbulent 
diffusion  of  microbubbles  by  inertial  subrange  turbulence  (Kolmogorov  turbulence). 

Hence  in  this  section  we  propose  that  di''.  wavenumber  spectrum  for  sound-speed 
fluctuations  in  the  bubble  layer  can  be  described  by  the  Kolmogorov  spectrum,  which  is 
proportional  to  K'^  1/3  ,  where  K  is  the  wavenumber  for  a  Fourier  component  of  the  sound- 
speed  distribution.  In  the  next  section  (Section  III)  we  use  the  spectrum  to  make  predictions 
for  scattering  from  the  bubble  layer. 

A.  Analysis  of  Reverberation  Data 

Figure  1  shows  some  typical  results  obtained  from  deei>ocean  reverberation 
experiments  for  scattering  versus  frequency  at  a  fixed  grazing  angle.  (Chapman  and  Harris, 
1962;  Chapman  and  Scott  1964;  Percy,  1970;  Brown  et  al.,  1966).  To  use  the  data  to 
estimate  the  horizontal  wavenumber  spectrum  P(qh)  as  a  function  of  qh.  we  convert  the 
measured  backscatter  strength,  BS  (in  decibels),  to  a  scattering  ci  )ss  section  per  unit  area, 

S ,  and  divide  by  the  geometric  factor  F  as  indicated  in  Eq.(33)  below 

P(<lh)=  ^  meu  /F  (33) 

where  the  measured  cross  section  ^nwa*  ‘s  related  to  the  backscatter  strength  BS  by  BS  = 
101ogio(4 

By  plotting  log,oF(qh)  versus  log,o  (qn)  from  the  measured  data,  we  observed  that, 
for  all  cases  examined,  the  wavenumber  spectrum  is  well  represented  by  a  power  law  of  the 
form  Aqh“P ,  where  A  is  a  strength  parameter  and  p  is  a  spectral  roll-off  exponent 


Numeiical  values  for  the  quantities  A  and  p  were  obtained  from  a  least-squares  fit  for  each 
data  set.  As  examples,  the  estimated  wavenumber  spectra  for  four  of  the  data  sets  in  Fig.  1 
are  shown  in  Fig.  2.  The  horizontal  axis  in  Fig.2  is  the  logarithm  ("base  10)  of  qh,  and  the 
vertical  axis  is  the  logarithm  (base  10)  of  the  estimated  wavenumber  spectrum,  P.  In 
determining  the  optimum  least-squares  fit,  the  e-folding  depth  was  varied  about  L  =1.5  to 
further  improve  the  fit.  Since  there  was  little  improvement,  in  all  four  cases  in  Fig.  2  we 
have  used  an  e-folding  depth  of  L=1.5  m,  a  value  that  is  consistent  with  oceanographic 
measurements  (Thorpe,  1982;  Farmer  and  Vagle  .1989;  Su  1992;  Canmill  and  Su,  1992; 

Su  et  al.,  1993). 

The  procedure  described  above  for  inferring  P(qh)  was  used  to  analyze 
reverberation  data  for  surface  backscatter  for  a  total  of  12  different  measurements,  for 
windspeeds  ranging  from  approximately  15  kts  to  over  30  kts.  In  Fig.  3  we  show  the 
infared  values  for  the  spectral  exponent  p  as  a  function  of  windspeed  for  the  12 
measurements  (Chapman  and  Hands,  1962;  Chapman  and  Scott  1964;  Percy,  1970;  Brown 
et  al„  1966;  Jin  et  al,  1989). 

In  a  given  experiment  the  strength  parameter  A  (and  hence  the  rms  sound-speed 
fluctuation  ac(0))  generally  increases  dramatically  and  sytematically  with  windspeed. 
However,  in  the  examination  of  the  12  measurements,  we  found  that,  between  different 
experiments,  the  correlation  of  the  strength  parameter  A  with  windspeed  was  poor.  Hence 
a  plot  of  the  strength  parameter  A  versus  windspeed  for  the  12  measurements  gives  little 
insight  and  consequently  is  not  presented. 

Although  the  strength  parameter  A  depends  strongly  on  windspeed,  the  spectral 
exponent  p  appears  to  have  little  windspeed  dependence.  The  mean  value  of  p  from  the 
12  least-squares  fits  is  p  =  3.86,  and  the  standard  deviation  is  .45. 

One  can  understand  the  origin  of  the  apparent  universal  value  p=  3.86  from  a 
simple  analysis  of  the  geometric  factor  F  and  an  inspection  of  plots  of  3  versus 
frequency  at  a  fixed  grazing  angle  (E.g.,  Fig.  1).  At  a  grazing  angle  of  2(P  ,  for  example. 


and  frequencies  above  a  few  kilohertz,  the  vertical  wavelength  is  on  the  order  of  the  e- 
folding  distance  ("thickness")  of  the  bubble  layer.  Consequently,  the  quantity  kj  in  the 
geometric  factor  is  much  greater  than  1,  so  that  the  frequency  dependence  in  the  geometric 
factor  comes  solely  from  the  k'J  factor  (See  Eq.(31)).  Thus  at  frequencies  above  a  few 
kilohertz,  the  geometric  factor  F  varies  as  frequency  to  the  fourth  power.  From  inspection 
of  plots  of  S  versus  frequency  at  a  fixed  grazing  angle,  we  can  see  that  <9  varies 
weakly  with  frequency  for  frequencies  above  a  few  kilohertz.  For  <9  to  depend  weakly  on 
frequency,  the  exponent  Ln  the  wavenumber  spectrum  P(qh)  niust  almost  cancel  the  fourth 
power  frequency  dependence  in  the  geometric  factor.  Therefore  we  can  directly  infer  from 
the  weak  frequency  dependence  of  <9  that  p  is  slightly  less  than  4  .  Since  <9  apparently 
always  has  weak  dependence  on  frequency  above  a  few  kilohertz,  the  data  direedy  imply  a 
nearly  universal  value  for  p  that  is  slighdy  less  than  4. 

The  consistency  of  the  value  of  p  for  the  12  measurements  suggested  to  us  that,  at 
the  scales  probed  by  the  backscattered  waves  (1/2  the  acoustic  wavelength),  the  distribution 
of  microbubbles  is  governed  by  a  model-independent  mechanism  such  as  small-scale 
turbulence.  In  the  discussion  below  we  present  evidence  for  the  turbulence  hypothesis. 

B.  Kolmogorov  Turbulence  and  the  Distribution  of  Microbubbles 

Langmuire  (1938),  in  his  original  studies  of  the  near-surface  circulation  that  now 
bears  his  name,  noted  that  "...although  the  motion  was  very  slow  it  was  very  turbulent." 

In  recent  years,  Thorpe  (1982)  has  made  comprehensive  observations  of  the  turbulent 
nature  of  Langmuire  circulation  and  emphasized  its  role  in  the  turbulent  diffusion  of 
microbubbles  (Thorpe  1982, 1984).  Very  recendy,  Osborn  et  al.  (1992)  simultaneously 
measured  both  the  turbulence  and  acoustic  scattering  associated  with  bubble  plumes  and 
concluded  that  turbulence  is  a  dominant  feature  of  the  plumes.  From  a  simplistic  viewpoint, 
the  bubble-cloud  records  from  upward-looking  high-frequency  sonars,  which  are  often 
remarkably  similar  to  clouds  in  the  atmosphere,  visually  suggest  the  turbulent  nature  of  the 


currents  that  transport  microbubbles  (Thorpe,  1982;  Farmer  and  Vagle,  1989;  Farmer, 
1992). 

Because,  at  small  scales,  turbulence  has  some  universal  characteristics,  we  seek 
here  to  establish  a  connection  between  the  empirical  value  of  p  =  3.86  and  small-scale 
turbulence.  We  first  present  some  basic  background  information  in  order  to  make  our 
reasoning  easier  to  follow. 

In  turbulent  fluids,  steady-state  small-scale  turbulence  is  governed  by  an  energy 
cascade  process  that  leads  to  model-independent  structure  at  small  length  scales  where  the 
fluid  has  "forgotten"  the  details  of  the  driving  mechanism  that  supplies  the  energy  (Landau 
and  Lifshitz,  1959;  Tatarski,  1961).  These  length  scales  are  commonly  known  as  the 
"inertial  subrange"  (Hinze,  1959;  Neumann  and  Pierson,  1966;  Grant  et  al.  1962;  Grant  et 
al.,  1968;  Embleton  and  Daigle,  1991 ).  Kolmogorov  has  shown  that  the  velocity  structure 
of  such  turbulence  has  a  three-dimensional  wavenumber  specmim  proportional  to 
s=  ,  regardless  of  the  intensity  of  the  turbulence  (Tatarski,  1961) .  Further, 

Obukhov  and  others  have  studied  the  the  concentration  distribution  of  "passive  additives" 
such  as  heat  (or,  in  our  case,  microbubbles)  that  are  mixed  by  Kolmogorov  turbulence 
(Tatarski,  1961;  Cbrrsin,  1951).  Obukhov  was  the  first  to  show  that,  given  sufficient  time, 
the  concentration  distribution  acquires  the  same  spectrum  as  the  turbulence  itself  (Tatarski, 
1961). 

In  the  ocean  mixed  layer,  Obukhov's  prediction  is  strongly  supported  by 
measurements  of  near-surface  temperature  fluctuations  (Whitmarsh,  et  al.,  1957;  Voorhis 
and  Perkins,  1966;  Grant,  et  al.,  1968)  which  clearly  show  the  characteristics  of  turbulent 
mixing  by  inertial  subrange  turbulence.  Whitmarsh,  et  al.,  for  example,  have  shown  that 
the  "structure  function"  for  the  temperature  fluctuations  <  [T(r2)  -  T(ri)]2>  closely 
follows  Kolmogorov’s  well  -  known  "2/3  law"  which  states  that  for  inertial  subrange 
turbulence  the  structure  function  varies  as  ,  where  p  =  -  ri|  is  the  separation 

between  two  measurement  points.  In  wavenumber  space  the  2/3  law  corresponds  to  a 


three-dimensional  wavenumber  spectrum  that  varies  as  (Tatarski,  1961).  Thus  four 
pieces  of  information  —  direct  observations  with  high-frequency  sonars, 
Kolmogorov/Obukhov  theory,  temperature  fluctuation  measurements,  and  the  acoustic 
analysis  presented  here  -  have  led  us  to  hypothesoe  that,  except  for  the  monotonic  decay 
with  depth  due  to  dissolution,  the  microbubble  distribution,  like  the  distributions  for 
turbulence  and  temperature,  can  be  represented  by  the  spectrum  associated  with 
inertial  sulnange  turbulence. 

The  inverse  power  law  K'^  1/^  proposed  by  Kolmogorov  is  not  valid  at  all  scales 
but  only  in  the  inertial  subrange,  2n/Lo  "<  K  «  2jt/lo  ,  v/here  L©  is  the  scale  of  the 
largest  eddies  ("outer  scale")  and  1^  is  the  scale  for  the  smallest  eddies  ("inner"  scale)  where 
energy  is  dissipated  by  viscosity.  The  complete  turbulence  spectrum  is  often  approximated 
with  an  analytic  expression  due  to  Von  Karman  (Clifford,  1978).  The  normalized  Von 
Karman  spectrum  varies  as  in  the  inertial  subrange  and  by  construction  has  a 

three-dimensional  integral  of  unity  in  wavenumber  space.  To  make  backscatter  predictions 
from  first  principles  without  a  direct  measurement  of  the  strength  parameter  A,  one  must 
use  some  expression  such  as  the  Von  Karman  spectn?  hat  can  be  normalized  in 
wavenumber  space.  Although  often  used  in  the  atmospheric  community,  the  Von  Karman 
spectrum  is  only  one  of  a  many  possible  choices  for  a  convenient  integrable  expression 
and  has  no  deep  physical  significance. 

For  our  purposes,  it  is  convenient  to  generalize  the  spectral  exponent  in  the  Von 
Karman  spectrum  and  set  the  inner  scale  to  zero.  With  a  general  spectral  exponent  p  and 
with  Iq  =  0,  the  normalized  Von  Karman  spectrum  is 

C(K)  =  N(P)  (8n3/2  Ko-3)(1  +  K2/Ko2)-P/2  (34) 

where  Ko  =  2r/Lo  ,  and  N(P)  =  F (P/2)/  F 0/2-3/2),  is  the  ratio  of  two  gamma 
functions.  In  the  inertial  subrange  we  have  K/Ko  »  1.  Hence  we  have 


<i(K)  =  N(p)  (8n3/2  Ko'^)(K/Ko)-P.  K/Ko  »  1  (35) 


Langmuire  observed  that,  in  Lake  George,  the  downwelling  currents  were 
consistently  confined  to  the  mixed  layer  (Langmuire,  1938).  Therefore,  in  the  ocean,  when 
a  mixed  layer  is  present,  we  expect  Lq  to  be  on  the  order  of  the  mixed  layer  depth,  or 
roughly  50  m  to  100  m.  This  estimate  is  consistent  with  existing  oceanographic 
measurements  of  horizontal  bubble  layer  structure  (Commander,  1992;  Henyey  et  al., 

1992).  To  get  some  idea  of  the  size  of  <i(K)  for  such  outer  scale  values,  we  can  consider 
the  special  case  p  =  4,  for  which  <i(K)  can  be  simply  evaluated.  For  p  =  4  and  L©  =  100 
m,  we  have 

<!:(K)  =  (16r2  /  Lo)K-^ 

(36) 

Any  value  of  p  that  is  close  to  4  gives  a  similar  numerical  coefficient  for  K'P. 

The  complete  wavenumber  spectrum  for  sound-speed  fluctuations  is  Oc^(O)  times 
the  normalized  spectrum  ^(K).  Hence,  using  the  result  in  Eq.(36),  we  have  for  the 
strength  parameter  A  (the  coefficient  of  K“P  inferred  from  reverberation  data)  the  value  A 
«  1 .6<yc^(0).  From  Fig.  2,  it  can  be  seen  that  for  windspeeds  in  the  range  of  20  kts  to  30 
kts,  the  inferred  values  of  logioA  (i.e.,  the  y^intercept)  are  in  the  range  3.5  to  5.2.  We  can 
therefore  infer  Oc(0)  to  be  in  the  range  of  44  m/s  to  3 15  m/s,  which  corresponds  to  rms 
void  firaction  fluctuations  of  order  10*^  to  10*5. 

The  results  for  ac(0)  inferred  from  backscatter  data  are  consistent  with  the  in-situ 
measurements  of  Su  (1992)  and  Cartmill  and  Su  (1992).  Using  acoustic  measurements  of 
bubble  sizes,  they  infer  values  for  <yc(0)  for  a  30  kt  wind  that  are  in  the  range  150  m/s  to 
200  m/s,  or  approximately  ten  percent  of  the  speed  of  sound  in  bubble-free  water. 
Preliminary  direct  travel-time  measurements  of  sound  speed  recently  reported  by  Lamarrc 


and  Melville  (1993)  for  moderate  sea  states  (winds  less  than  20  kt)  also  suggest  near¬ 
surface  sound  speed  fluctuations  on  the  order  of  ten  percent  Farmer  and  Vagle  (1989),  in 
contrast,  have  estimated  near-surface  average  sound  speed  reductions  of  about  1  per  cent 
If  we  assume  that  the  fluctuations  in  sound  speed  are  on  the  same  order  of  magnitude  as  the 
average  sound-speed  reduction,  then  our  inferred  value  for  Cc(0)  and  the  measured  values 
of  Su,  Cartmill  and  Su,  and  Lamaire  and  MelviUe  are  all  significantly  larger  than  the  value 
inferred  from  the  data  of  Farmer  and  Vagle.  In  the  backscatter  predictions  in  this  article  we 
use  data  for  05(2)  provided  by  Su  (1992)  since  they  are  the  only  in-situ  measurements 
presently  available  for  high  sea  states. 

It  should  be  noted  that  although  the  Von  Karman  spectrum  used  here  is  isotropic  at 
all  scales,  isotropy  is  a  plausible  approximation  only  in  the  inertial  subrange.  Even  in  the 
inertial  subrange,  measurements  of  temperature  fluctuations  show  a  clear  departure  from 
isotropy  (Whitmarsh,  et  al.,  1957).  Fortunately,  it  is  easy  to  generalize  the  scattering 
model  to  account  for  anisotropy  in  the  normaliz»l  wavenumber  spectrum  0(K).  For  low 
grazing  angle  backscatter,  the  primary  anisotropy,  the  vertical/horizontal  anisotropy,  affects 
only  the  normalization  of  the  spectrum  and  not  the  spectral  exponent  p  .  To  account  for 
vertical/horizontal  anisotropy ,  we  simply  multiply  the  normalized  wavenumber  spectrum 
for  isotropic  turbulence  by  a  factor  Lv/^h .  where  Ly  and  Lh  are,  respectively,  the  outer 
scales  associated  with  vertical  and  horizontal  turbulence  structure.  Without  oceanographic 
measurements  of  Ly  and  Lh ,  however,  any  attempt  to  theoretically  account  for  anisotropy 
would  be  pure  conjecture.  Hence,  in  this  article,  we  normalize  the  spectrum  under  the 
assumption  of  isotropy  (Ly  =  Lh=  Lo)  realizing  that  such  an  assumption  is  crude  at  best. 
Ideally  we  would  like  to  have  a  direct  measurement  of  the  spectrum,  but,  currently,  no 
measurements  exist 


C.  RMS  Sound-Speed  Fluctuations  Versus  Depth 

To  complete  the  statistical  description  of  the  bubble  layer,  we  need  an  estimate  of 
the  rms  sound-speed  fluctuation  as  a  function  of  depth.  As  discussed  in  Section  I,  we  use 
an  exponential  fitting  function  of  the  form  Oc(2)  =  Oc(0)e~*/L ,  an  approximation  that  is 
consistent  with  observed  bubble  distributions  (Thorpe,  1982  ;  Farmer  and  Vagle,  1989; 

Su,  1992;  Su  et  al.,  1993).  For  the  backscatter  predictions  in  a  30-kl  wind  presented  in  the 
next  section,  the  function  Cc(z)  was  fitted  to  the  rms  sound-speed  measurements  of  Su 
(1992)  from  a  measurement  site  in  the  Northeast  Pacific  in  winter  conditions.  A  least- 
squares  fit  to  the  measured  values  of  versus  depth  gives  an  rms  sound-speed  fluctuation 
at  the  surface  of  Oc(0)  =  207  m/*  and  an  e-folding  distance  of  L=i  .95.  In  the  next  section, 
the  Htted  values  for  <Jc(0)  and  L,  together  with  the  above  result  for  ^:(K)  are  used  in 
Eqs.(30)-(32)  to  confute  the  backscatter  cross  section. 


ni.  COMPARISON  OF  BACKSCATTER  PREDICTIONS  WITH  EXPERIMENT 

In  the  previous  section  we  proposed  a  wavenumber  spectrum  for  use  in  a  stochastic 
scattering  model.  In  this  section  we  use  the  proposed  spectrum  and  the  stochastic  scattering 
model  to  make  predictions  for  backscatter  from  the  oceanic  bubble  layer.  We  compare  the 
predictions  with  the  experimental  results  repoiial  by  Ogden  and  Erskine  (1992)  for  the 
Critical  Sea  Tests  (CST).  The  purpose  of  the  comparison  is  to  show  that  with 
occanognq)hically  constrained  parameters  (no  "adjustable"  parameters)  the  stochastic  model 
can  give  a  good  account  of  the  measured  backscatter  as  function  of  both  frequency  and 
angle. 

We  consider  two  sets  of  calculations.  The  first  set  uses  the  empirically  determined 
value  of  ^  3.86  and  the  second  uses  the  Kolmogorov  value  for  the  inertial  subrange,  p  = 
11/3  =  3.67.  In  both  calculations  we  use  the  normalized  Von  Karman  spectrum  with  an 
outer  scale  of  Lo  =  106  m,  which  was  the  average  mixed  layer  depth  for  a  30  kt  wind 
during  the  CST-7  experiment  (Farmer,  1992).  As  discussed  in  the  previous  section,  the 
rms  sound-speed  fluctuation  profile  was  fitted  to  measured  data  of  Su  (1992)  for  a  30  kt 
wind. 

Figure  4  compares  the  predicted  scattering  strength-vs-grazing  angle  at  several 
firequencies  with  the  Ogden-Erskine  empirical  fits  to  CST  data.  Since  at  very  low 
frequencies,  surface  scattering  can  be  the  dominant  scattering  mechanism,  we  have  for 
completeness  included  two  theoretical  curves,  one  that  includes  both  bubble  layer  scattering 
and  rough  surface  scattering(solid  line),  and  one  with  only  bubble  layer  scattering(dashed 
line).  Tlie  rough  surface  contribution  is  from  Thorsos  (1990)  and  is  added  incoherently  to 
the  bubble-layer  contribution.  Figure  5  shows  the  backscatter  as  a  function  of  the 
frequency  for  sev^al  different  fixed  grazing  angles.  Since  the  solid  and  dashed  lines  in 
Figs.  4  and  5  overlay  except  at  vay  low  frequencies,  we  can  see  that  the  main  scattering 
contribution  is  from  the  bubble  layer. 


Figures  6  and  7  are  the  same  as  Figs  4  and  5  except  the  Kolmogorov  v<due  of 
P=1 1/3=3 .67  was  used  instead  of  p=3.86.  With  ^3.67  the  agreement  with  the  CST  data 
is  not  quite  as  good  as  with  p=3.86,  but  nevertheless  is  satisfactory  and  shows  that 
Kolmogorov's  theory  for  inertial  subrange  turbulence  is  consistent  with  the  observed 
backscatter. 


IV.  SUMMARY  AND  CXDNCLUSIONS 

A  stochastic  scattering  model  based  on  the  wavenumber  spectrum  for  sound-speed 
fluctuations  in  the  oceanic  bubble  layer  has  been  derived.  The  model  was  used  to  analyze 
data  from  12  different  measurements  of  surface  backscatter.  From  the  analysis  it  was  found 
that  the  effective  horizontal  wavenumber  spectrum  for  sound-speed  fluctuations  in  the 
bubble  layer  is  an  inverse  power  law  of  the  form  P(K)=AK‘P,  where  A  is  a  strength 
parameter,  K  is  the  horizontal  wavenumber,  and  the  mean  value  of  the  spectral  roll-off 
exponent  is  p  =  3.86  ±  .45  .  The  universal  character  of  p  suggests  that  the  small-scale 
structure  of  the  microbubble  distribution  (scales  <  5  -  10  m)  is  governed  by  turbulent 
mixing  due  to  turbulence  in  the  inertial  subrange.  Moreover,  using  a  spectrum  based  on  the 
notion  of  inertial  subrange  turbulence  and  observed  oceanographic  parameters,  good 
agreement  was  obtained  with  measurements  from  the  Critical  Sea  Test  surface  backscatter 
experiments. 

Although  the  stochastic  scattering  model ,  with  no  adjustable  parameters,  gives 
good  agreement  with  surface  backscatter  data,  there  nevertheless  is  a  need  for  further 
oceanographic  measurements  to  establish  the  physical  validity  of  the  model.  Since  the 
backscatter  depends  on  the  square  of  the  rms  sound-speed  fluctuation,  the  most  important 
experiment  is  a  direct  sound-speed  measurement  to  confirm  the  measurements  of  Su  (1992) 
where  Oc(z)  was  inferred  indirectly  from  acoustic  estimates  of  the  bubble  size  distribution. 
The  next  most  important  measurement  is  to  detCTmine  whether  or  not  an  inertial  subrange 
actually  exists  in  the  spectrum  of  the  sound-speed  fluctuations.  Establishing  that  p  is 
exactly  1 1/3  is  not  as  important  as  determining  whether  the  bubble  distribution  is 
systematically  controlled  by  small-scale  turbulence.  In  fact,  since  the  value  1 1/3  applies  to 
fully  developed  isotropic  turbulence,  it  would  not  be  surprising  if  the  effective  value  for  6 
in  the  bubble  layer  were  not  exactly  1 1/3.  Finally,  an  important  and  hopefully  possible  task 
would  be  to  relate  the  observed  wavenumber  spectrum  for  sound-speed  fluctuations  in  the 
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bubble  layer  to  environmental  parameters  such  as  windspeed,  air/sea  temperatures,  and  sea 
state.  To  accomplish  such  a  task  would  require  that  the  wavenumber  spectrum  and  the 
associated  environmental  parameters  be  measured  under  a  variety  of  conditions  so  that 
predictive  models  could  be  tested. 

We  have  presented  evidence  that  scattering  from  microbubbles  in  the  near-surface 
bubble  layer  is  as  strong  as  any  other  surface  scattering  mechanism,  but  we  have  not 
shown  that  it  is  the  sole  mechanism.  If,  for  example,  scattering  from  dense,  near-surface 
bubble  plumes  under  breaking  waves  were  equally  strong,  the  cross-section  would  be 
raised  by  3  dB,  an  amount  that  is  within  the  experimental  error  in  the  measurements.  Thus, 
until  more  definitive  measurements  are  made,  the  stochastic  model  presented  here  must  be 
viewed  as  one  of  a  number  of  mechanisms  that  can  contribute  significantly  to  surface 
backscatter. 
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HGURE  CAPTIONS 


Fig.  1.  Backscatter  strength  versus  frequency  for  a  grazing  angle  of  2(P.  Note  the  common 
weak  dependence  on  frequency  at  frequencies  above  a  few  kilohertz. 

Fig.  2.  Horizonal  wavenumber  spectra  for  sound-speed  fluctuations  in  the  bubble  layer. 

The  spectra  were  inferred  from  the  backscatter  versus  frequency  curves  in  Fig.  1.  The  form 
of  the  spectra  is  P(qh)  =  A  qh'P  where  A  is  a  strength  parameter  and  p  is  a  spectral  roll-off 
exponent.  The  quantity  qh  is  2kh ,  where  kh  is  the  horizontal  wavenumber. 

Fig.  3.  The  spectral  roll-off  exponent  ^  versus  windspeed  for  12  different  backscatter 
measurements.  There  appears  to  be  litde  dependence  on  windspeed.  The  approximately 
constant  value  of  p  =  3.86,  is  a  result  of  the  weak  frequency  dependence  seen  in  Fig.  1. 

The  nearly  universal  value  of  p  suggests  that  the  wavenumber  spectrum  of  sound-speed 
fluctuations  in  the  bubble  layer  is  governed  by  a  model-independent  mechanism  such  as 
inertial  subrange  turbulence  (Kolmogorov  turbulence)  for  which  p  =  1 1/3  =  3.67  for  fully 
developed  isotropic  turbulence. 

Fig.  4.  Scattering  strength  versus  grazing  angle  at  several  fixed  frequencies  using  the 
empirical  value  of  p=3.86.  The  dashed  line  is  backscatter  just  from  the  bubble  layer.  The 
solid  line  is  the  incoherent  sum  of  bubble  layer  scattering  and  rough  surface  scattering.  The 
diamonds  are  the  Ogden-Erskine  empirical  fit  to  the  Critical  Sea  Test  data. 

Fig.  5.  Same  as  Fig.  4  except  that  the  grazing  angles  are  fixed  and  the  frequency  is  varied. 

Fig.  6.  Same  as  Fig.  4  except  the  Kolmogorov  value  of  p  =  1 1/3  =  3.67  is  used. 

Fig.  7.  Same  as  Fig.  4  except  the  grazing  angles  are  fixed  and  the  frequency  is  varied.  Also 
the  Kolmogorov  value  of  p  =  1 1/3  =  3.67  is  used. 
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Transducer  Development 

Principal  Investigator:  M.  A.  Breazeale 


Research  Accomplished  in  1992: 

Most  of  the  year  was  spent  on  transducer  development  and  in  completing  the  Ph.D. 
Dissertation  of  Dehua  Huang;  however,  some  of  the  time  was  spent  in  completing  research 
supported  in  previous  years  under  Navy  Programs.  These  will  be  included  in  this  Final  Technical 
Report. 

Ph.D.  Dissertation  of  Dehua  Huang.  The  Ph.D.  Dissertation  "Gaussian  finite  element 
method  for  description  of  underwater  sound  diffiraction,"  was  completed  in  August  1992.  The 
dissertation  describes  a  new  method  for  solving  diffraction  problems.  The  method  is  based  on  the 
use  of  Gaussian  diffraction  theory.  The  Rayleigh  integral  is  used  to  prove  the  core  of  Gaussian 
theory:  the  diffraction  field  of  a  Gaussian  radiator  also  is  described  by  a  Gaussian  function.  The 
parabolic  approximation  used  by  previous  authors  is  not  necessary  to  this  proof.  Comparison  of 
the  Gaussian  beam  expansion  and  Fourier  series  expansion  reveals  that  the  Gaussian  expansion  is 
a  more  general  and  more  powerful  technique. 

The  method  combines  the  Gaussian  beam  superposition  technique  [Wen  and  Breazeale,  J. 
Acoust.  Soc.  Am.  83,  1752-1756  (1988)]  and  the  Finite  element  solution  to  the  parabolic 
equation  [Huang,  J.  Acoust.  Soc.  Am.  84,  1405-1413  (1988)].  Computer  modeling  shows  that 
the  new  method  is  capable  of  solving  for  the  sound  field  even  in  an  inhomogeneous  medium, 
whether  the  source  is  a  Gaussian  source  or  a  distributed  source.  It  can  be  used  for  horizontally 
layered  interfaces  or  irregular  interfaces.  Calculated  results  are  compared  with  experimental 
results  by  use  of  a  recently  designed  an  improved  Gaussian  transducer  in  a  laboratory  water  tank. 
In  addition,  the  power  of  the  Gaussian  Finite  element  method  is  demonstrated  by  comparing 
numerical  results  with  experimental  results  from  use  of  a  piston  transducer  in  a  water  tank.  A 
publication  detailing  these  results  is  in  preparation. 


Patent  Application.  A  patent  application  has  been  filed  to  cover  the  advances  made  in  the 
fataicarion  of  a  Gaussian  transducer.  The  patent  in  the  names  of  Dehua  Huang  and  M.  A. 
Breazeale  has  been  filed  by  Airmar  Technology  Corporation,  current  employer  of  Dehua  Huang, 
on  December  8,  1992.  It  is  anticipated  that  a  number  of  licenses  will  be  issued  for  various 
applications  of  the  Gaussian  principle  to  different  situations. 

Elacidc  Nonlinearities  in  Single  Crystal  Gallium  Arsenide.  Experimental  investigation  of 
single  crystal  gallium  arsenide  nonlinearity  has  been  completed  and  the  results  have  been 
published.  Gallium  arsenide  is  a  technologically  important  crystal  whose  nonlinearity  had  not 
been  investigated.  The  hannonic  generation  technique  is  ideally  suited  for  measurement  of  a  set 
of  third  order  elastic  constants  of  this  material,  and  the  variation  of  ultrasound  velocity  with 
applied  hydrostatic  pressure  completes  the  set  of  six  third  order  elastic  constants.  By  using  a 
modification  of  the  Keating  theory  for  cubic  lattice  solids  we  were  able  to  obtain  a  complete  set  of 
third  order  elastic  constants  between  77K  and  300K. 

Nonlinear  Techniques  for  Nondestructive  Evaluation  of  Composites.  A  new 
investigation  has  been  opened  up  with  the  use  of  nonlinear  techniques  to  investigate  composites 
and  heat  damage  in  composites.  Since  very  little  scientific  infonnation  is  available  about  the 
physical  behavior  of  composites,  it  is  highly  desirable  that  new  techniques,  such  as  the  harmonic 
generation  technique,  be  applied  to  collection  of  data  on  composites.  It  also  desirable  to  use  more 
standard  techniques  to  collect  data  on  composites.  Tentative  measurements  have  been  made  of 
the  behavior  of  composites.  Tentative  measurements  have  been  made  of  the  behavior  of 
composites.  We  find  that  the  velocity  in  the  basal  plane  is  approximately  one-third  of  that  in  the 
axial  direction.  Thus,  composites  are  highly  anisotropic  the  behavior  of  composites;  however,  it  is 
clear  that  primitive  cells  in  composites  are  not  hexagonal.  Thus,  much  is  to  be  learned  about  the 
linear  and  the  nonlinear  physical  properties  of  composites. 
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ELASTIC  nonlinearities  IN  SINGLE  CRYSTAL 
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Third  order  elastic  <TOE)  constants  arcinecessary  (o  a 
fundamental  description  of  many  physical  jiropertics  of  solids 
resulting  from  lattice  nonlinearity.  The  properties  in  quesdon 
include  thermal  expansion,  heat  conduction,  tempenture 
dependence  of  the  specific  heat,  temperature^  4od;  pressure 
dependence  .of  the  elutic  constants,  piffeteoce^pveen  the 
adiabatic  and  isothermal  elastic  constant^'  phonon  yis^ity, 
thermal  attenuation  of  acoustic  wave^  shock  .Reformation  of 
materials,  acoustically  induced. static  strpss  j^^strain,.etc. 
Previously'-^  we  have  shown  that  measureifay'it  of  ultrasonic 
harmonic  generation  gives  informatioo  about  cotdbinations  of  TOC 
constants  from  helium  temperatures  to  at  least  3^*  iC  Since  no 
single  technique  gives  all  six  TOE  constants  of  cubic^crystals  it 
was  necessary  to  introduce  a  second  technique  in  older  to  plot  all 
six  of  the  TOE  constants  of  the  semiconductor  Gallium  Arsenide 
(GaAs)  over  a  wide  tenqieraure  range.  The  purpose  of  this  report 
is  to  describe  the  technique  used  because  it  makes  optimum  use  of 
theory  and  experiment  to  arrive  at  data  which  would  not  be 
available  otherwise.  The  technique  shruld  be  applicable  without 
modification  to  the  evaluation  of  all  TOE  constants  of  all 
zincblende  structure  compounds.  The  data  on  GaAs  verify  its 
validity. 


The  propagation  of  a  finite  amplitude  ultrasonic  wave 
along  the  thrw  (xrinci-'al  directions  in  a  cubic  lattice  is  described  by 

^a2u  ,8u  a^u  ... 


Assuming  an  initially  sinusoidal  wave  allows  one  to  write 
a  solution  in  the  fotm 


U  =  A 1  sin  (ka  -  ox)  +  A2  cos  2(ka  -  o«)  +  ... 


A2  =  -l3  +  |^)Ajk2x.  (3) 

Defining  the  nonlinearity  parameter 

P  =  -13  +  1^)  (4) 

allows  one  to  write  the  nonlinearity  parameters  in  terms  of 
measured  quantities 

P  =  •  <5) 

A  k2  X 
I 


where  k  =  —  is  the  propagation  constant  and  x  is  the  sample 

length.  In  terms  of  TOE  constants  the  nonlinearity  parameters 
along  the  principal  directions  are 

Pl00  =  '[3  +  ^1 


*110  =  (3 


2(C|  1  +  Ci2  +  2C44) 


±I2Ci66. 

+  2C44) 


(6  a,b.c) 


Pm  =13+ 


3(C|  I  +  2Ci2  +  4C44) 


All 


By  measuring  the  amplitudes  A2,  Ai,  the  frequent^  (from  which  k 
can  be  calculated),  and  the  sample  length  x  in  the  principal 
directions  one  can  evaluate  the  noolinearify  parameters.  The 
results  for  GaAs  are  shown  in  Fig.  1  plotted  between  77‘  R  and 
room  temperature. 


TEMPERATURE  (K) 

Hg.  1.  Temperature  dependence  of  nonlinearity  parameten  of 
GaAs. 

To  evaluate  the  individual  TOE  constants  from  these  data, 
additional  information  is  needed.  The  Keating^  theory  is  a  three- 
parameter  theory  from  which  one  can  calculate  idl  six  TOE 
constants.  Thus,  in  principle  these  data  plus  the  Keating  theory 
would  be  sufficient  to  isolate  ail  six  TOE  constants  of  GaAs.  We 
chose  to  provide  additional  cxperimenial  data,  however. 


The  most  accurate  evaluation  of  TOE  constants  of  cubic 
crystals  at  room  temperature  appears  to  come  from  the 
combination  of  harmonic  generation  data  with  those  taken  by  use 
of  pressure  variation  of  ultrasonic  wave  velocity.*  A  pressure 
bomb  was  used  to  take  pressure  variation  data  at  room  temperature 
with  the  same  GaAs  samples.  A  plot  of  the  normally  Cr^uency 
(essentially  sound  velocity)  as  a  function  of  pressure  is  ^ven  in 
Fig.  2.  The  curves  all  are  straight  lines  except  for  the  longitudinal 
wave  in  the  [111]  direction  Oabclled  [111],[111)  in  Fig.  2).  The 
slopes  of  these  curves  can  be  used  to  calculate  combinations  of 
TOE  constants  at  room  temperature. 
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Fig.  2.  Variation  of  normalized  frequency  (ultrasonic  wave 
velocity)  with  pressure  in  GaAs. 


Results 

The  results  of  the  two  sets  of  measurements  can  be 
combined  to  isolate  all  six  TOE  constants  at  room  lempcrature. 
The  results  have  been  evaluated  for  GaAs  and  compared  with 
room  temperature  values  of  other  researchers.  Having  these 
values  one  now  is  able  to  use  the  Keating  model  along  with  the 
harmonic  generation  data  to  calculate  the  temperature  dependu^ 
of  each  TOE  constant.  Results  of  the  values  of  all  six  TOE 
constants  of  GaAs  between  room  temperature  and  Uquid  nitrogen 
temperature  are  given  in  Figs.  3  and  4  in  which  the  curves  are  least 
squi^res  fits  of  the  data  with  a  fifth  order  polynomial.  Most  of  the 


TOE  constants  are  linear  functions  of  temperature;  however,  both 
Cn2  and  Ci23  exhibit  remarkable  temperature  variations  and 
emphasize  the  importance  of  beina  able  to  measure  these 
fundamental  parameters  c^er  a  range  of  tempoatuies. 


Fig.  3.  TempcraturcdependenceofTOEconstantsCiii.Ciae 


tndCas*. 


T  n  M  P  H  R  A  T  U  R  E  (  K  ) 

Fig.  4.  Temperature  dependence  of  Cm.  Cm  and  C144 
CaiK-hv  Rrl«t.nr^^ 

Although  the  Cauchy  relations  usually  are  evaluated  for  the 
«ocond  order  elastic  constants  (from  the  sound  velocity),  they  also 
can  be  evaluated  for  the  TOE  constants  as  well  If  all  lattice 
mteractions  arc  of  the  central  force  type,  and  each  atom  ts  a  center 
of  symmetry,  then  all  elastic  constants  should  obey  the  Cauchy 
nelatior.s.  For  the  TC3E  constants  the  Cauchy  relations  arc 

Cll2  =  C!55  (7) 

C123  =  Ci«4  =C456  (8) 

If  in  addition  the  shon  range  repulsive  interaction  is  the  dorrunant 
one,  there  are  the  so-called  "strong  '  Cauchy  xlations 


We  have  examined  the  strong  (Zauchy  relations  over  the  available 
temperature  range  by  comparing  the  measured  quantities 

Cl  12  4Ci55  ttnd  C123  +  6Ci44  +  SCasfi  with  |Cin  and  find 

that  there  is  a  tendency  for  better  agreement  with  the  TOE  constant 
Cauchy  relations  as  OK  is  approached.  Since  both  the  Keating 
model  and  the  Ciuchy  relations  strictly  should  apply  only  at  OK, 
diis  may  indicated  that  the  lattice  interaction  in  GaAs  is  basically  of 
the  ceiiirnl  type,  and  the  deviations  from  the  Qiuchy  relations  are 
caused  by  thermal  effects.  An  analogous  observation  has  been 
made  for  germanium  and  silicon.^ 

Summary 

All  six  TOE  constants  of  GaAs  have  been  evaluated 
between  room  temperature  and  liquid  nitrogen  temperatures.  Two 
of  the  TOE  constants,  Cn2  and  C123,  exhibit  considerable 
variation,  (at  least  a  I'sctor  of  10  for  C123).  The  other  four  arc 
almost  linear  functions  o''  temperature. 

The  technique  wc  have  used  can  be  applied  directly  to 
evaluation  of  the  TOP  constants  of  all  zincblende  structure 
compounds.  Other  structures  or  samples  in  which  interstitials  or 
dislocations  are  prominent  require  further  analysis. 
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INTRODUCTION 

When  I  was  asked  to  review  the  progress  in  nonlinear  acoustics  1  realized  that 
considerable  progress  recently  has  been  made  even  in  a  relatively  restricted  field,  the 
nonlinear  acoustics  of  solids.  For  example,  there  recently  was  a  correlation  of  the 
nonlinearity  parameter  in  zincblende  structure  crystals  with  interatomic  potential  functions.’ 
Later  a  generalization^  showed  that  in  essentially  all  symmetries  of  crystalline  solids  one 
finds  that  the  nonlinearity  parameter  depends  exclusively  on  the  Bom-Mayer  hardness 
parameter,  meaning  that  in  all  crystalline  solids  the  magnitude  of  the  nonlinearity  parameter 
is  largely  determined  by  the  shape  of  the  interatomic  potential  function.  Furthermore,  with 
crystdline  solids  having  a  larger  nonlinearity  parameter  one  finds  an  increasingly  imporpint 
zero-point  energy  contribution.^  The  nonlinearity  parameter  also  has  been  correlated  with 
the  Anderson-Griineisen  parameter  in  solids^  and  temperature  dependence  of  the  third 
order  elastic  constants  of  NaCl-structure  alkalidc  halide  crystals.^  Of  more  technological 
interest  is  the  fact  that  the  nonlinearity  parameter  has  been  conelated  with  the  ultimate  yield 
strength  of  isotropic  solids^  and  with  hardness  in  steels.^  Nowadays  the  magnitude  of  the 
nonlinearity  parameters  of  such  inhonx^cneous  substances  as  muscle  tissue”,  PZT 
ceramics^  and  high  Tc  superconductors^such  as  YBa2Cu307-6  ,  arc  being  measured. 
Results  on  some  of  these  materials  form  the  substance  of  this  report. 

With  ihc  wide-ranging  investigations  to  be  reported,  one  tends  to  give  the 
impression  that  there  has  been  randomness  in  the  development  of  nonlinear  acoustics,  that 
it  "just  grew"  as  my  title  implies.  The  truth  is  that  it  did.  Nonlinear  acoustics  just  grew. 

As  early  as  1660,  Hooke  discovered  that  we  could  get  reasonable  results  just  by  ignoring 
nonlinearity.  Subsequently,  more  subtle  observations  required  a  nonlinear  theory.  In  this 
case  there  still  was  the  desire  to  "keep  things  simple"  and  use  only  a  linear  theory.  Whether 
or  not  the  "linearization"  worktd  depended  on  the  detail  desired  in  the  agreement  between 
experiment  and  theory.  At  that  time  the  understanding  of  nonlinearity  was  not  sufficient  to 
allow  one  to  predict  where  and  how  it  would  be  encountered.  This  is  the  reason  for  the 
lack  of  systematic  development.  Early  on.  scientists  did  not  understand  that  the  nonlinear 
theory  is  the  more  fundamental  one.  Having  become  enlightened,  some  scientists  now  are 
developing  a  systematic  approach  to  nonlinearity,  but  many  of  the  developments  have  been 
random,  analogous  to  the  way  Topsy  in  Uncle  Tom's  Cabin  "just  grew."  Harriet  Beecher 
Stowe  wrote  Uncle  Tom's  Cabin^in  Stowe  House  near  Bowdoin  College,  so  it  is 
appropriate  that  1  quote  directly  from  her  dcscTiption  of  Topsy  .so  we  can  have  a  mental 
picture  of  how  Topsy  grew,  and  by  analogy,  how  nonlinear  acoustics  grew. 
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One  morning,  while  Miss  Ophelia  was  busy  in  some  of 
her  domestic  cares,  St.  Clare’s  voice  was  heard,  calling  her 
at  the  foot  of  the  stairs. 

"Come  down  here,  cousin;  I've  something  to  show 
you.” 

"What  is  it?",  said  Miss  Ophelia,  coming  down,  with 
her  sewing  in  hand. 

"I’ve  made  a  purchase  for  your  department,  —  see 
here,"  said  St.  Clare;  and,  with  the  word,  he  pulled  along  a 
little  Negro  girl,  about  eight  or  nine  years  of  age. 

Sitting  down  before  her,  she  began  to  question  her. 

"How  old  arc  you,  Topsy?" 

"Dunno,  Missis,"  said  the  image,  with  a  grin  that 
showed  all  her  teeth. 

"Don’t  know  how  old  you  are?  Did  n’t  anybody  ever 
tell  you?  Who  was  your  mother?" 

"Never  had  none!"  said  the  child,  with  another  grin. 

"Never  had  any  mother?  What  do  you  mean?  Where 
were  you  bom?" 

"Never  was  bom!"  persisted  Topsy,  with  another  grin, 
that  looked  so  goblin-like,  that,  if  Miss  Ophelia  had  been  at 
all  nervous,  she  might  have  fancied  that  she  had  got  hold  of 
some  sooty  gnome  from  the  land  of  Diablerie;  but  Miss 
Ophelia  was  not  nervous,  but  plain  and  business-like,  and 
she  said,  with  some  sternness,  — 

"You  mus  n’t  answer  me  in  that  way,  child;  I  ‘m  not 
playing  with  you.  Tell  me  where  you  were  born,  and  who 
your  father  and  mother  were." 

"Never  was  bom  reiterated  the  creature,  more 
emphatically;  ‘never  had  no  father  nor  mother,  nor  nothin’. 

I  was  raised  by  a  speculator,  with  lots  of  others.  Old  Aunt 
Sue  used  to  take  car  on  us." 

The  child  was  evidently  sincere;  and  Jane,  breaking  into 
a  short  laugh,  said,  -- 

"Laws,  Missis,  there  ‘s  heaps  of  ’em.  Speculators  buys 
’em  up  cheap,  when  they’s  little,  and  gets  'em  raised  for 
market." 

"How  long  have  you  lived  with  your  master  or 
mistress?" 

"Dunno,  Missis." 

"Is  it  a  year,  or  more,  or  less?" 

"Dunno,  Missis." 

"Laws,  Missis,  those  low  Negros,  -  they  can’t  tell; 
they  don’t  know  anything  about  time,"  said  Jane;  "they  don’t 
know  what  a  year  is;  they  don’t  know  their  own  ages." 

"Have  you  ever  heard  anything  about  God,  Topsy?" 

The  child  looked  bewildered,  but  grinned  as  usual. 

"Do  you  know  who  made  you?" 

"Nobody,  as  I  knows  on,"  said  the  child,  with  a  short 
laugh. 

The  idea  appeared  to  amuse  here  considerably;  for  her 
eyes  twinkled,  and  she  added,  — 

"I  spect  I  grow’d.  Don’t  think  nobody  never  made 

me." 

This  is  my  impression  of  the  subject  of  nonlinear  acoustics.  1  have  yet  to  find 
somebody  who  claims  to  be  the  father  of  the  subject,  and  1  doubt  that  1  will  meet  the 
mother.  The  subject  “just  grow’d. ’’  Almost  twenty  years  ago  1  made  .some  observations 
for  the  International  Journal  of  Nondestructive  Testing.  *0  The  conclusion  written  at  that 
time  bears  repeating,  for  I  think  it  will  set  the  tone  for  this  entire  session. 
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Nonlinear  beliavior  of  solids  in  which  finite  amplitude 
ultrasonic  waves  propagate  can  be  demonsuated.  It  can  also 
be  predicted  from  an  extension  of  elasticity  tlicory.  It  is 
hoped  that  some  of  the  phenomena  observed  can  serve  as  tlic 
basis  of  new,  exceptionally  sensitive  nondestructive  testing 
techniques. 

Since  that  time  major  changes  have  taken  place  in  the  subject.  Nondestructive 
Testing  has  become  nondestructive  Evaluation.  In  addition  correlation  is  being  observed 
between  nonlinearity  and  work  hardening,  for  example,  or  nonlinearity  and  hardness  of 
steels.  Such  efforts  arc  of  fundamental  importance  to  nondestructive  evaluation  and  to 
development  of  nondestructive  evaluation  techniques.  There  is  another  aspect  of  the 
development  of  a  science  that  cannot  be  ignored.  That  is  the  basic  knowledge  of  nonlinear 
material  interactions  and  the  mathematical  description  of  them.  This  is  the  reason  we 
introduced  third  order  clastic  constants  in  the  early  sixties.  At  the  time  they  were 
introduced  we  didn't  have  the  foggiest  notion  about  their  meaning.  We  didn't  even  know 
that  third  order  clastic  constants  are  almost  always  negative  and  have  a  magnitude  of  the 
order  of  ten  times  the  magnitude  of  the  second  order  constants.  We  still  arc  trying  to 
understand  them,  but  we  know  more  about  them  than  we  did  then.  Nowadays  we  know 
enough  to  make  measurement  of  the  nonlinear  behavior  of  a  solid  and  try  to  correlate  the 
measurement  with  intrinsic  properties.  This  can  succeed,  but  it  will  succeed  only  to  the 
extent  that  we  thoroughly  understand  that  some  approximations  and  assumptions  are 
involved,  and  recognize  them  at  the  time  we  are  interpreting  data.  For  example,  there  arc 
such  things  as  intrinsic  nonlinearities  and  then  there  arc  nonlinearitics  arising  from  the 
characteristics  of  individual  samples.  The  intrinsic  nonlinearities  arise  from  interatomic 
forces  in  the  crystalline  lattice,  and  can  correctly  be  described  by  third  order  elastic 
constants.  The  other  nonlinearitics  can  arise  from  strains,  dislocations  or  other 
imperfections  in  the  sample.  Sometimes  our  measurement  techniques  will  separate  the  two, 
and  sometimes  not.  It  is  up  to  us  to  make  a  correct  interpretation  of  the  data.  Only  then  can 
we  develop  genuinely  dependable  techniques  for  NDE  based  on  sample  nonlinearity. 

Some  of  the  early  experiments  are  illustrative  of  what  we  have  to  be  aware  of  if  we 
are  to  make  sense  of  nonlinearity  measurements,  so  I  would  like  to  turn  back  history  for 
more  than  a  quarter-century  and  present  measurements  made  at  that  time,  to  show  how 
much  we  have  progressed  in  measurement -even  of  samples  we  measured  initially,  and  to 
end  with  recent,  very  provocative  measurements  on  piezoelectric  ceramics,  and  even  high 
Tc  superconductors. 

ARCHAEOLOGY  OF  NONLINEARITY 

If  we  start  digging  we  can  discover  that  one  participant  in  the  early  prehistory  of 
nonlinear  acoustics  was  my  colleague  Don  Thompson.  He  and  I  will  reminisce  while  the 
others  can  marvel  that  we  had  any  inkling  of  the  meaning  of  nonlinearity  so  far  back  in 
history.  In  order  to  facilitate  discussion,  it  is  necessary  to  write  down  the  basic  form  of  the 
nonlinear  equation  and  of  its  solution.  Near  the  end  of  the  discussion  1  want  to  point  out 
some  implications  about  the  form  of  the  nonlinear  equation  that  comes  from  our  analysis  of 
PZT. 

THEORY 

To  correctly  derive  the  equation  desenbing  propagation  of  an  ultrasonic  wave  in  a 
nonlinear  crystalline  lattice  we  begin  with  the  definition  of  tlie  elastic  potential  energy  in 
terms  of  strains  by  writing  an  expansion  in  the  fomi: 


•K^)  -  L  Qjkl  Rij  Hkl  -I,  X  Cij|,;lmn  RijRklRmn  d 

ijkl  ’  ijklmn 

The  Cijkl  are  the  clastic  constants  that  appear  in  the  linear  approximation.  Ihe  Cikjlmn  are 
the  third  order  clastic  constants.  They  arc  the  .set  of  coefficients  that  make  the  use  of  the 
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first  nonlinear  terms  meaningful.  At  this  point  one  can  introduce  tlic  effect  of  any 
additional  strain,  such  as  piezoelectricity,  etc.,  by  adding  appropriate  terms.  Once  we 
proceed  beyond  this  point,  however,  the  correct  introduction  of  additional  strains  becomes 
increasingly  difficult.  One  favorite  technique  for  forcing  nonlinear  acoustics  to  grow  (like 
Topsy,  to  be  sure),  is  to  include  only  the  elastic  terms  given  in  1  and  to  proceed  with 
the  derivation  of  the  wave  equation  in  which  one  defines  a  nonlinearity  parameter.  The 
next  step  is  to  see  whether  this  wave  equation  agrees  with  experiment  in  various  solids  by 
determining  whether  the  nonlinearity  parameter  has  in  it  effects  that  the  theory  does  not 
account  for.  If  it  docs,  one  should  begin  again  and  account  for  the  newly  discovered 
strains  in  Eq.  1  and  rcdcrivc  everything.  This  is  the  point  at  which  many  physicists  behave 
as  speculators.  They  try  to  market  their  results  rather  than  to  correctly  rederive  them.  To  a 
certain  extent,  the  authors  of  this  treatise  will  behave  as  a  speculator,  as  it  develops.  To 
compensate,  there  is  an  implied,  promise  to  complete  the  job  at  some  further  date.  Thus,  the 
speculator  has  integrity  after  all. 

There  are  several  ways  to  proceed  once  Eq.  1  has  been  written.  The  most  direct  is 
to  define  the  Lagrangian  function  and  use  Lagrange's  equation  to  derive  the  nonlinear  wave 
equation.  By  using  the  appropriate  form  of  Lagrange's  equations  and  specializing  to  a 
specific  orientation  of  the  coordinates  with  respect  to  the  ultrasonic  wave  propagation 
direction  one  can  write  the  nonlinear  wave  equation  in  the  form 


P  Xi 


—  i 

9aik=i  Unite/ 


(2) 


which  shows  exactly  how  the  strain  energy  with  the  elastic  constants  enters  the  nonlinear 
wave  equation.  If  the  strain  energy  correctly  accounts  for  action  of  all  of  the  forces  in  the 
sample,  then  the  derived  nonlinear  equation  correctly  describes  the  propagation  of  an 
ultrasonic  wave  in  the  sample.  Aye,  there's  the  rub.  1  am  not  sure  that  always  the 
averaging  effect  of  all  of  the  molecules  propagating  the  wave  is  adequately  described  by  the 
intrinsic  third  order  elastic  constants.  But  we  have  no  alternative  at  the  moment  but  to 
assume  that  it  does  and  proceed  with  our  analysis. 

The  thing  that  pleased  us  back  in  the  dark  ages,  and  the  thing  that  made  possible 
much  of  the  progress  in  nonlinear  acoustics,  was  the  fact  that  when  one  specializes  Eq.  2  to 
pure  mode  directions  one  finds  that  the  equations  are  separable  and  that  the  longitudinal 
wave  is  described  by  an  equation  of  the  form 


^  ^  ,a^ 

Po  — —  =  K2 - +  (3K2  +  K3) - 

at2  aa2  aa  aa^ 


(3) 


in  all  three  principal  directions  in  a  cubic  lattice,  and  in  an  isotropic  medium  as  well.  We 
wrote  the  equation  in  this  fomi  to  emphasize  the  fact  that  it  applies  in  all  principal 
directions,  then  defined  K2  and  K3  as  shown  in  Table  1  so  one  could  see  the  role  played  by 
the  third  order  elastic  constants.  The  definitions  of  the  nonlinearity  parameter 

P  (4) 

emphasizes  the  fact  that  it  is  the  ratio  of  the  coefficients  of  the  nonlinear  term  to  the  linear 
term  in  the  nonlinear  wave  equation.  TTiis  nonlinearity  parameter,  with  the  definitions 
given  in  Table  I,  only  accounts  for  clastic  nonlinearitics.  To  account  for  other  effects,  the 
expressions  for  the  nonlinearity  parameters  would  be  more  complicated. 

If  one  derived  the  nonlinear  wave  equation  in  a  slightly  different  way  he  could 
define  the  analogous  nonlinearity  parameter  for  fluids.  It  was  interesting  to  discover  that 
the  nonlinearity  parameter  for  solids  is  of  the  same  order  of  magnitude  as  the  nonlinearity 
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Tabic  1  K2 and  K3  for (100),  (110), and  (111)  Directions  in  Cubic  Crystals 


Direction 

K2 

K3 

(100) 

ClI 

Cm 

(110) 

Cn  +  C12  +  2C44 

Cm  +  3Cii2  +  12Cir)6 

2 

4 

(111) 

Cl  1  +  2ci2  +  4c44 

Cii  1  +  6C112  +  12Ci44  +  24Ci66 

3 

9 

2Ci23  +  I6C456 

9 

Table  II  Comparison  of  structure,  bonding  and  acoustic  nonlinearity  parameters  along  the 
[100]  direction  of  cubic  crystals. 


STRUCTURE 

BONDING 

RANGE  OF  P 

2Lincblendc 

Covalent 

1.8  -  3.2 

Flourite 

Ionic 

3.4  -  4.6 

FCC 

Metallic 

b 

< 

b 

FCC 

(inert  gas) 

Van  der  Waals 

5.8  -  7.0 

BCC 

Metallic 

5.0  -  8.8 

NaCl 

Ionic 

13.5  -15.4 

parameter  for  fluids.  This  can  be  confirmed  by  comparing  Tables  II  and  III.  Table  II  gives 
nonlinearity  parameters  for  a  number  of  solids.  Table  111  gives  B/A  and  P  =  B/A  +  2  for 
fluids.  For  solids  the  nonlinearity  parameters  range  from  2  to  15.  For  fluids  the  range  is 
from  6  to  14,  approximately  the  same.  This  range  can  serve  as  reference  for  some  data  I 
will  present  later. 

The  solution  c."  Cq.  3  allows  one  to  show  how  nonlinearity  parameters  can  be 
measured.  A  perturbation  solution  takes  the  form 

U  =  Ai  sin(ka  -  (ot)  -  P  A^  k^a  cos  2(ka  -  cot)  +  ...  (5) 

This  solution  shows  that  an  iiiitially  sinusoidal  ultra.sonic  wave  in  a  solid  will  produce  a 
second  harmonic  whose  amplitude  is  proportional  to  the  nonlinearity  parameter.  Thus,  if 
we  can  measure  the  amplitude  of  the  fundamental  and  tlic  second  harmonic  after  the 
ultrasonic  wave  has  propagated  through  the  sample,  we  caji  determine  the  nonlinearity 
parameter 


P  = 


_ 

A^  k^a 


(6) 


since  the  propagation  constant  k  =  — ,  where  c  is  the  wave  velocity.  The  .sample  length 
is  a.  Measurement  of  P  as  a  function  of  temperature  can  be  quite  informative. 
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Tabic  III  Values  of  B/A  and  P  for  fluids  at  atmospheric  pressure 


Liquid 

Temperature 

(*C) 

B/A 

P  =  B/A  +  2 

Water-Distilled 

0 

4.16 

6.16 

20 

4.96 

6.96 

40 

5.38 

7.38 

60 

5.67 

7.67 

80 

5.96 

7.96 

100 

6.11 

8.11 

Acetone 

30 

9.44 

11.44 

Benzene 

30 

9.03 

11.03 

Benzyl  Alcohol 

30 

10.19 

12.19 

CCI4 

30 

11.54 

13.54 

TEMPERATURE  (K) 

Figure  1.  Temperature  variation  of  the  nonlinearity  parameters 
in  the  principal  directions  of  NaCl. 


COMPARISON  OF  THE  NONLINEARITY  PARAMETERS  OF  NaCl,  PZT  and 
YBa2Cu307^ 

Recently  we  measured  the  nonlinear  behavior  of  NaCl.  The  nonlinearity  parameters 
for  the  principal  directions  in  NaCl  covered  a  wide  range  of  values,  but  each  did  not  change 
much  with  temperature  between  room  temperature  and  liquid  nitrogen  temperature. 
Measured  nonlinearity  parameter  values  in  the  three  principal  directions  in  NaCl  are  given 
in  F’3  1.  Different  samples  were  used  for  each  direction;  however,  measured  values 
differed  only  in  the  1 1 1 1  ]  direction.  The  Oiigin  of  the  discrepancy  is  not  known,  even 
though  it  originally  was  thought  to  be  the  effect  of  OH'  ions.  It  is  clear  that  the  nonlinearity 
parameter  in  the  (1 1 1  ]  direction  is  extremely  sensitive  to  sample  difference,  however. 
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The  sample  of  NaCl  is  a  single  crystal,  so  there  is  variation  in  the  magnitude  of  the 
nonlinearity  parameter  for  different  orientations.  When  one  has  a  sample  such  as  lead 
zirconate  titanate  (yZV)  one  has  a  ceramic  which  doesn’t  necessarily  exhibit  crystalline 
properties.  We  recently  measured  the  nonlinearity  parameters  of  two  types  of  PZT  both  in 
the  polarized  state  and  in  the  unpolarizcd  state.  The  nonlinearity  paiameter  was  measured 
along  the  direction  of  polarization.  Results  for  the  two  samples  of  PZT  are  given  in  Figs. 
2,  3, 4,  and  5  as  a  function  of  temperature.  In  the  cubic  crystal  NaCl,  the  magnitude  of  the 
nonlinearity  parameter  is  14  or  less.  Well  below  the  Curie  temperature  Tc,  the  nonlinearity 


Tc 


Temperature  (°Q 

Figure  2.  Temperature  dependence  of  the  nonlinearity  parameter  in 
K1 -Unpolarizcd  sample  of  PZT. 


Tc 


Temperature  {°C) 

Figure  3.  Temperature  dependence  of  the  nonlinearity  parameter  in 
K 1  -Polarized  sample  of  PZT. 
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Tc 


Temperature  (^C) 

Figure  4.  Temperature  dependence  of  the  nonlinearity  parameter  in 
Sl-Unpolarizcd  sample  of  PZT. 


Tc 


Figure  5.  Tcnipcraturc  dependence  of  the  nonlinearity  parameter  in 
SI -Polarized  sample  of  PZT. 


parameters  in  the  samples  of  PZT  are  in  this  range.  However,  as  the  Curie  temperature  is 
approached  the  magnitude  of  the  nonlinearity  parameter  becomes  anomalously  high,  as 
high  as  1500  in  one  case.  The  origin  of  this  anomalously  high  nonlinearity  parameter  is  not 
known  at  the  moment.  We  assume  that  the  large  nonlinearity  parameter  is  produced  by 
effects  other  than  elastic  ones.  Thus,  it  becomes  necessary  for  this  speculator  to  modify 
Eq.  1  by  including  other  considerations.  Once  this  is  done,  it  will  be  ncccs.sary  to  rcdcrive 
the  nonlinear  wave  equation.  Only  then  will  one  be  certain  whether  the  anomalously  large 
magnitudes  of  the  nonlinearity  parameters  of  PZT  really  have  meaning,  or  whether  a 
completely  new  theory  is  necessary  to  explain  this  new  dcvclopnKnt  in  the  field  of 
nonlinear  acoustics. 
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Figure  6.  Temperature  variation  of  P  in  the  high  Tc  superconductor 
YBa2Cu307-«(Tc'=94-K). 


The  final  example  is  the  behavior  of  tlie  nonlinearity  parameter  of  the  high  Tc 
superconductor  YBa2Cu307.s.  We  measured  the  second  harmonic  as  a  function  of 
temperature  through  the  transition  temperature  Tc  between  the  normal  state  and  the 
superconducting  state.  When  we  interpreted  the  results  in  terms  of  the  nonlinearity 
parameter  p,  we  got  the  results  shown  in  Fig.  6.  At  room  temperature  (30"C)  the 
nonlinearity  parameter  had  a  large  value,  approximately  14.  At  lower  temperatures  the 
attenuation  became  noticeable,  so  we  corrected  for  attenuation.  The  result  is  that  both  in  the 
corrected  data  and  in  the  uncorrected  data  the  nonlinearity  parameter  appears  to  vanish  at  the 
transition  temperature  Tc,  which  is  just  below  lOO'K.  This  anomalous  effect  in  the 
behavior  of  the  ceramic  YBa2Cu307^  is  just  opposite  to  that  of  PZT  at  the  Curie 
Temperature.  Considerable  theoretical  speculation  now  is  desirable  to  explain  the  behavior 
of  the  nonlinearity  parameters  of  ceramics  near  transition  temperatures. 
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V  Temperature  dependence  of  elastic  nonlinearities  in  single-crystal 
gallium  arsenide 

D.  Joharapurkar  •’  D.  Gerlich,^”  and  M.  A.  Breazeale'’ 

National  Center  for  Physical  Acoustics.  University  of  Mississippi,  University,  Mississippi  38677 
(Received  27  August  1991;  accepted  for  publication  22  May  1992) 

The  six  third-order  clastic  moduli  (TOEM)  of  single-crystal  gallium  arsenide  were  determined 
by  a  combination  of  measurements  of  ultrasonic  second-harmonic  generation,  and  pressure 
dependence  of  the  second-order  elastic  moduli,  at  room  temperature.  In  the  temperature  range 
77-3(X)  K,  the  nonlinearity  parameter  for  the  propagation  directions  [100],  [1 10],  and  [1 1 1]  was 
measured.  Utilizing  the  Keating  model,  these  data  were  used  in  evaluating  all  six  TOEM  as  a 
function  of  temperature.  The  TOEM  Cm,  and  C456  turn  out  to  be  nearly  constant  in  the 
above  temperature  range.  The  Cauchy  relations  seem  to  be  obeyed  somewhat  better  as  0  K  is 
approached.  The  measured  values  of  the  TOEM  have  been  employed  in  calculating  a 
Mumaghan  equation  of  state,  which  predicts  a  somewhat  higher  volume  change  than  the 
measured  one.  The  elastic  Griineisen  constants  deduced  from  the  TOEM  are  in  reasonable 
agreement  with  the  thermal  ones  in  the  high-ten.perature  limit. 


1.  INTRODUCTION 

The  nonlinear  elastic  properties  of  diamond  and  zinc- 
blende  structure  materials  has  been  the  subject  of  numer¬ 
ous  studies.'  In  most  of  these  studies  only  the  room- 
temperature  properties  of  the  materials  have  been 
investigated,  and  only  relatively  few  have  been  measured  as 
function  of  temperature.'"^  Several  investigators  have  de¬ 
termined  the  third-order  elastic  moduli  (TOEM)  of  GaAs 
in  the  past,^  but  all  this  work  has  been  restricted  to  room 
temperature.  Since  the  elastic  nonlinearities  arc  often  cor¬ 
related  with  various  other  anharmonic  properties,  knowing 
the  temperature  dependence  of  the  elastic  nonlinearities  is 
quite  vital  in  many  instances.  Due  to  the  lack  of  more 
detailed  information,  researchers  sometimes  have  assumed 
that  the  elastic  nonlineariiies  in  solids,  specifically  the 
TOEM,  are  temperature  independent.’  In  order  to  clarify 
this  question  somewhat  further,  the  present  study  was  un¬ 
dertaken  to  measure  the  TOEM  of  single-crystal  GaAs  as 
a  function  of  temperature.  We  have  measured  all  six 
TOEM  at  room  temperature  by  combining  results  of 
second-harmonic  generation  (SHG)'  and  hydrostatic  pres¬ 
sure  dependence  of  the  second-order  elastic  moduli 
(SOEM).*  In  the  temperature  range  77-300  K,  the  three 
linear  combinations  of  the  TOEM  available  from  SHG 
have  been  determined  experimentally.  Utilizing  the  Keat¬ 
ing  model,’  and  its  extension  to  the  zinc-blende  structure,'® 
all  six  TOEM  have  been  determined.  This  is  possible  be¬ 
cause  the  validity  of  the  extended  Keating  model  at  room 
temperature  can  be  verified  by  experimental  data.  Thus, 
this  investigation  presents  all  six  TOEM  of  GaAs  as  a 
function  of  temperature  in  the  range  77-300  K. 
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The  Keating  model’  has  been  spectacularly  successful 
in  describing  the  linear  and  anharmonic  elastic  properties 
of  the  diamond  structure  materials  germanium  and  silicon. 
Since  the  number  of  adjustable  parameters  in  the  model 
(two  for  the  second  order,  three  for  the  third  order)  is 
smaller  than  the  number  of  elastic  moduli  in  each  case 
(three  second  order  and  six  third  order),  the  applicability 
of  the  model  may  be  readily  tested  by  comparing  with 
experimental  data.  From  measurement  of  SHG  in  the  three 
longitudinal  pure  propagation  modes  ([100],  [110],  and 
[111])  ill  the  temptrature  range  4-300  K,  and  utilizing  the 
Keating  model,  all  six  TOEM  as  a  function  of  temperature 
could  be  determined."  It  should  therefore  be  of  interest  to 
examine  whether  the  Keating  model’  and  its  extension'® 
are  applicable  to  the  zinc-blende  structure  as  well,  and 
whether  all  six  TOEM  as  a  function  of  temperature  may  be 
evaluated  in  the  same  way.  These  qucstioiu  arc  examined 
in  the  present  work. 


II.  EXPERIMENT 

Single-crystal  boules  of  GaAs,  n-type.  Si  doped,  having 
a  room-temperature  resistivity  of  1.1  X  10“'  fl  cm,  carrier 
concentration  of  1.2  X  10‘*  cm“',  carrier  mobility  of  1990 
cm’ V“'  s“',  were  purchased  from  Crystal  Specialties  In¬ 
ternational.  Two  right  parallelepipeds  were  cut  from  the 
above  boules,  the  faces  of  one  corresponding  to  crystalline 
planes  ( 100),  ( 1 10),  and  ( iTO),  the  other  one  correspond¬ 
ing  to  ( 1 1 1 ),  ( iTO),  and  (112)  planes.  Opposite  faces  were 
lapped  flat  and  parallel.  For  the  SHG  measurements,  lon¬ 
gitudinal  ultrasonic  waves,  fundamental  frequency  of  30 
MHz,  propagating  in  the  [100],  [1 10],  and  [111]  crystalline 
directions  were  utilized  For  the  hydrostatic  pressure  runs, 
longitudinal  and  shear  waves,  15  MHz  frequency,  propa¬ 
gating  in  the  [110]  and  [111]  directions  were  used.  The 
waves  were  generated  with  quartz  crystalline  transducers, 
X  and  Y  cut,  bonded  to  the  specimens  for  the  room- 
temperature  measurements  with  benzephenone,  while  for 
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the  li)w-tcmpcraturc  measurements  the  backing  of  scotch 
tape  was  used  as  a  bonding  agent. 

The  room-temperature  ultrasonic  velocity  as  well  as  its 
change  under  hydrostatic  pressure  were  determined  by  the 
gated  double-pulse  superposition  method/^  which  pro¬ 
vided  a  resolution  of  10“  .  The  hydrostatic  pressure  was 
generated  by  means  of  a  small  manual  hydraulic  pump, 
using  light  weight  oil  as  the  pressurizing  fluid,  generating  a 
pressure  range  of  0-0.2  GPa.  The  sample  was  placed  inside 
a  stainless-steel  pressure  vessel  which  was  connected  to  the 
pump  and  a  Bourdon  gauge,  which  served  for  measuring 
the  pressure  in  the  system.  These  measurements  yield  the 
values  of  the  three  SOEM,  C||,  c,2,  and  c^^,  as  well  as  the 
three  linear  combinations  of  the  TOEM,‘^ 

C'i  =  Qii  +  2Cii2. 


Q  =  Q44  +  2Ci55,  (1) 

i 


In  the  SHG  measurements,  one  determines  the  ampli¬ 
tudes  of  the  fundamental  and  second-harmonic  wave, 
and  Alt  and  these  are  directly  related  to  the  nonlinearity 
parameter  P,  which  is  given  by 


/3=8 


3^2+ Ki 


(2) 


where  k  is  the  wave  vector,  a  the  length  of  the  sample. 
and  Kj  are  the  second-  and  third-order  coupling  constants. 
These  constants  are  linear  combinations  of  the  SOEM  and 
TOEM,  respectively,  e.g.,  in  the  cubic  system  (Laue  group 
Cl),  for  the  three  pure  propagation  modes,  [100],  [110], 
and  [111],  they  are  given  by' 

/r2|ioo]=c,„ 


^2!  110]  —  (cji  -)-Ci2-(-2c44)/2, 


(3) 


^2!  H 1 1  —  (<^ii  +  2ci2-f-4c44)/3. 


and 


/rj[ioo]=c„„ 


1 10]  —  (Cjii  -f-3C|  i2-(-  12Ci55)/4,  (4) 


^3[  1 1  1  ]  —  (Qll  +  6Qi2  +  2Ci23-(-  12Ci44-f  24Ci55 


-I- 16C45«)/9. 

As  can  be  seen  from  Eqs.  (1)  and  (4),  measuring  the 
pressure  derivatives  of  the  SOEM  together  with  SHG  mea¬ 
surements  for  the  [100],  [ilO],  and  [111]  direction  enables 
one  to  determine  of  all  six  TOEM  with  their  propagated 
errors.  The  six  TOEM  are  as  follows: 

c,„=/:j(ioo], 

c, 23=^3(1001-03, 

c,  /r,(iooi 

c,,2-  2  y— . 

^  AjIllU)  A3IIOU)  c, 


„  „  2a:3[iio]  /:3[ioo]  c, 

144  2  3  12  4’ 


^  9A:3[111]  3C,  3C2  C3 

16  “  16  ~  4 

The  corresponding  errors  are  as  follows: 

|ac,„|  =  |aa:,[ioo]|, 

|AC,23|<|A/:3[100]|-f-|AC3l, 

,  |AC,|  .  |AAr,[100]| 

I  AC,,2|  <  -  H - j - , 


,  ia/:3[iio]|  1aa:3[100]|  |ac, 
I^QssI  <  ^  1-  + 


|AC,44|<|AC2H- 
.  I  AC,  I 


3  '  24  ■  8 

2|aa:3[iio]|  |Ai:3[ioo] 


12 


and 


9|AA:3[111]|  3|AC,|  3|AQ1  IAC3I 


8 


In  the  present  investigati^'n,  at  room  temperature  the 
values  of  for  the  three  directions  were  evaluated  by 
measuring  the  absolute  values  of  Ai  and  Ai-  In  the  low- 
temperature  run,  a  relative  measurement  of  P  was  carried 
out.  This  was  effected  by  keeping  both  values  of  A^  and  Ai 
constant  during  the  temperature  run.  A ,  was  kept  constant 
by  varying  the  intensity  of  the  sound  energy  delivered  to 
the  transducer,  while  A2  was  kept  constant  by  changing  the 
bias  voltage  on  the  capacitive  transducer.  As  can  be 
shown,'  the  following  relation  is  obtained  for  relative  val¬ 
ues  of  P: 


p{T)  _VtiT„)Co{T„)K2{T) 

P(Tk)  V,{T)Co{T)K2{T^)  -  ^  ’ 

•Here,  T  is  the  temperature,  Tg  a  reference  temperature 
(normally  room  temperature),  and  Cp  is  the  capacity  of 
the  capacitive  detector  utilized  in  the  determination  of  A^ 
and  A 2-'*  The  temperature  dependence  of  the  SOEM,  re¬ 
quired  for  deducing  the  values  of  from  P,  was  taken 
from  the  work  of  Cottam  and  Saunders.”  We  measured 
sound  velocities  in  the  three  principal  directions  and  found 
that  our  values  agree  within  experimental  uncertainty  with 
those  of  Cottam  and  Saunders. 


!ll.  RESULTS 

For  the  pressure  dependence  of  the  SOEM,  five  inde¬ 
pendent  propagation  modes  were  available:  one  longitudi¬ 
nal  and  two  shear  modes  in  the  [110]  direction,  and  a 
lonvitiidina!  and  shear  mode  in  the  fill]  direction.  The 
expenmoiiLal  setup  \ieius  a  lesouam  irequcuuv  y,  as  a 
function  of  pressure  P.  In  Fig.  1  we  present  the  changes  in 


j.  I.  Normalized  resonance  frequency  as  a  function  of  pressure  for 
ious  propagation  modes  at  room  temperature:  •  (1 10],  (1 10];  +(1 1 1], 
I];  A  [MO],  (001];  ■  [III],  [OlT];  4  [110],  [iTO]. 


:  normalized  resonant  frequency  ///o  as  a  function  of  P, 
lere  /o  is  the  resonant  frequency  at  P—Q.  The  dots  are 
r  experimentally  measured  data  points,  while  the  straight 
es  are  a  least-squares  fit  to  these  points.  The  slopes  of 
:se  lines,  d{  f^f{^/BP,  are  directly  related  to  the  pres- 
re  derivatives  of  the  natural  modulus,  where  pg  is 
:  equilibrium  density  and  W  is  the  natural  velocity," 

(6) 

:re,  the  subscript  0  denotes  the  equilibrium  state,  while 
:  prime  denotes  the  pressure  derivative,  (pof^^)o  the 
'EM  for  the  specified  direction  and  is  denoted  by  lu.  The 
ues  of  (po"  )o'  *s  deduced  from  the  slopes  of  Fig.  1, 
lether  with  their  associated  errors,  are  shown  in  Table  I. 
e  error  bars  for  the  individual  data  points  in  Fig.  1  are 
nificantly  smaller  than  the  statistical  scatter.  Hence  the 
tistical  scatter  of  the  least-squares  fit  to  the  straight  line 
ised  as  the  probable  error.  In  Table  11,  the  values  of  the 
sar  combination  of  the  TOEM,  C^,  Cj,  and  Cj,  together 
h  their  errors,  are  shown,  as  well  as  the  analogous  quan- 
es  derived  from  the  work  of  MeSkimin  and  Andreatch.’ 
;ure  2  presents  plots  of  A 2  vs  a]  for  the  three  SHG 
ipagation  directions  [100],  [110],  and  [1 1 1]  at  room  tem- 
ature.  The  dots  are  the  data  points,  while  the  straight 
s  are  least-squares  fits  to  the  points.  For  these  data  at 
5t  seven  measurements  of  each  sample  showed  that  the 
roducibility  is  better  than  6%.  The  slopes  of  these  lines, 
t/BA^,  are  presented  in  Table  III.  From  the  slopes  A",  is 
:ulated  for  the  principal  directions  and  the  results  are 


3LE  I.  Values  of  (Pd**^^)o  1°''  vanous  propagation  modes 


le 

(110] 

(1 10]  shear 

1 1 10]  shear 

(llij 

1111] 

long 

(iTO]  pol 

(001)  po! 

long 

shear 

4  9g 

*0.050 

-0.090 

*0.016 

0  750 

±0  038 

4  90 

±0  29 

0  1 30 
±0.049 

TABLE  II.  Values  of  the  linear  combinations  of  the  TOEM,  C,,  Cj,  and 
Cj  (units  are  GPa). 


c, 

c, 

Present  work 

-1401*1 

-512*8 

-538*8 

MeSkimin  and  Andreatch* 

-1396 

-536 

-565 

“Reference  5 


shown  in  Table  IV,  together  with  similar  data  derived  from 
the  work  of  MeSkimin  and  Andreatch.’  As  can  be  seen,  the 
agreement  of  both  sets  of  data  with  the  results  of  Me¬ 
Skimin  and  Andreatch  is  very  good.  From  the  data  shown 
in  Tables  II  and  IV,  the  six  room-temperature  TOEM  may 
now  be  determined,  and  the  results  are  shown  in  Table  V. 
Also  shown  in  the  same  table  are  analogous  results  ob¬ 
tained  by  other  investigators.^  As  can  be  seen,  the  agree¬ 
ment  among  the  different  results  for  all  six  room- 
temperature  TOEM  is  very  reasonable. 

In  Figs.  3  and  4,  the  nonlinearity  parameter  /3  and  A3 
are  shown  as  a  function  of  temperature  in  the  range  77-300 
K.  As  can  be  seen,  both  0  and  A3  vary  less  than  25%  over 
this  temperature  range,  although  a  temperature  depen¬ 
dence  is  obvious.  For  every  sample  the  experiment  was 
repeated  at  least  three  times  for  temperature  runs  to  pro¬ 
duce  the  data  plotted. 

IV.  DISCUSSION 
A.  Keating  model 

The  Keating  model’  has  been  very  successful  in  de¬ 
scribing  the  clastic  properties  of  diamond  structure  mate¬ 
rials.  The  theory  is  very  plausible,  physically  intuitive,  re¬ 
quires  a  minimal  number  of  adjustable  parameters,  and  is 
cast  in  a  rotationally  invariant  form.  It  expresses  the  qua¬ 
dratic  part  of  the  strain  energy  in  terms  of  two  second- 
order  force  constants:  a  bond-stretch  constant  a,  and  an 
angle-bend  consunt  P.  The  cubic  part  is  expressed  in  terms 
of  three  third-order  force  constants;  a  bond-stretch  con¬ 
stant  y,  an  angle-bend  constant  5,  and  a  mixed  bond- 


FIG  2  A,  as  a  funclion  of  A\  from  ihc  rtxmi-icmpcraiure  SHG  mcasurc- 
mc  nis 
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TA^LE  III.  V»Iues  of  (he  slopes  for  differenc  propagation  di¬ 

rections  (units  are  lO’ m”'). 


TABLE  IV.  The  room-temperature  values  o( P,  A'j.  and  Aj  deduced  from 
SHG  measurements  (A'j  and  A'3  are  in  units  of  GPa) 


[100] 

[110] 

[111] 

1.591:*:  0.025 

1.865  ±0.016 

0.958  ±0.035 

Stretch  angle-bend  constant  e.  The  model  has  been  ex¬ 
tended  to  the  zinc-blende  structure,'®  and  in  this  case  the 
three  SOEM  and  the  six  TOEM  are  given  by  the  following 
expressions; 


f^ii 


a-h3ff 

a 


4.053 


a-P 

c,2  = - -5.538  -T-. 


4a/3 


Cyn=r-^  +  ^€+'^1-201{ZW/a*), 

Cn2  =  r-5-l-e+1.531(ZV/a"). 

Ci23  =  y  4-  35  -  3e  -(-  24.663  ( ZV  V ) . 

c,«=y(i-^')4-fi(i+?')+e(i+?)(3^-i) 

a-P  ,  ,  ZW 

-f ( 24. 663  -  3  3 . 526^  -  0. 820^^ ) 


or 

C,55  =  y(l-^)^-6(l+f^)+e(l-ff)(3-^)+ - 

a 

1  -^V 

-I-  ( 1.531  -33. 526f-f  10.062^^)-^. 

Q56  =  y(l-5:)^-(- (24.663-50.288^ -f42.753f^ 
-19.203^^)-^, 


Propagation 

direction 

A': 

e 

A, 

A,‘ 

[100] 

1 18.4 

2.3 

-628  ±54 

-622 

[110] 

145  2 

5.62 

-1251±7 

-1252 

[1111 

154.2 

4.19 

-I108±23 

-  1 124 

‘Reference  5. 


Coulombic  long-range  interaction.  In  Table  VI  the  short- 
range  and  Coulombic  contributions  to  the  various  TOEM 
at  room  temperature  are  shown.  As  can  be  seen,  except  for 
C|44  and  Ci55  and  possibly  C123,  the  long-range  Coulombic 
contribution  is  completely  negligible. 

Since  there  are  only  three  adjustable  parameters  for  the 
six  TOEM  in  the  Keating  model,  from  the  three  measured 
values  of  ATj  one  may  determine  all  three  third-order  force 
constants  y,  5,  and  £,  and  thus  calculate  all  six  TOEM.  In 
Fig.  5,  the  values  of  y,  5,  and  f  as  a  function  of  temperature 
are  shown  as  calculated  from  the  data  of  Fig.  It  is  in¬ 
teresting  to  note  that  while  5  and  e  are  nearly  constant  over 
the  whole  temperature  range,  y  exhibits  a  considerable 
variation.  In  Figs.  6  and  7,  the  values  of  the  TOEM  as 
calculated  from  the  third-order  force  constants  6,  c,  and  y 
in  Fig.  5  are  shown.  The  curves  are  a  fifth-order  polyno¬ 
mial  fit  of  the  data.  The  TOEM  C||2,  C123,  and  C155  show 
considerable  variation  over  the  temperature  range  mea¬ 
sured.  The  other  three  TOEM  are  nearly  linear  functions 
of  temperature  in  the  range  77-300  K.  A  comparison  of  the 
values  of  the  TOEM  in  Table  V  and  those  read  from  Figs. 
6  and  7  is  a  good  test  of  the  validity  of  the  Keating  model 
at  room  temperature.  Such  a  comparison  shows  that  all  of 
the  TOEM  in  Figs.  6  and  7  except  C144  are  within  the 
errors  stated  in  Table  V.  Closer  examination  of  the  values 
for  Cm  given  by  the  diCerent  authors  listed  in  Table  V 
shows  a  wide  range  of  measured  values.  C144  appears  to  be 
the  most  uncertain  of  the  TOEM. 


( (a-^)/a]  -  10.058(ZV/a‘') 
((a+^)/a]-4.189(ZV/o*)  ’ 

Here,  a  is  the  lattice  constant  of  the  unit  cell,  Z  the  effec¬ 
tive  charge  number,'*  and  q  che  electronic  charge.  The 
terms  containing  Z^q^/a*  are  the  corrections  due  to  the 


B.  Cauchy  relations 

If  all  lattice  interactions  are  of  central  force  type,  and 
each  atom  is  a  center  of  symmetry,  the  elastic  moduli 
should  obey  the  Cauchy  relations.  For  the  SOEM  and 
TOEM,  these  relations  are  as  follows: 


TABLE  V  Room-tcmpcraiurc  values  of  the  TOEM  (units  arc  GPa) 


1  ! 

C,u 

^173 

c,„ 

^456 

Present  work 

-628  ±54 

-  387±  27 

-  90  ±  62 

24  1  ]  7 

-  269  ±  5 

-  44  ±  20 

MeSkimin  and 

Andreatch* 

-627 

.187 

-  ?7 

T 

269 

- 

Drabble  and 

Brammer'’ 

-  67^ 

4fX1 

4 

'■(5 

^20 

Abe  and  Imai' 

-620 

Vi 

!  4 

:k: 

44 

■Reference  5 
'Tteference  4 
'Reference  6 
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FIG.  3.  Nonlinearity  parameter  ^  as  a  function  of  temperature. 

=  ^112  =  ^155,  C^2i  =  C^u  =  C4i6■  (8) 

If,  in  addition  to  the  above  conditions,  the  short-range  re¬ 
pulsive  interaction  is  the  dominant  one,  there  are  addi¬ 
tional  relations  between  the  elastic  moduli,'^  the  so-called 
“strong  Cauchy  relations,”  which  are 

^11  =  ^12  =  2^4,,, 

C||i  =  2Cu2=2Cij5,  (9) 

Q23  =  ^144  =  Qji  =  0- 

From  Eqs.  (8)  and  (9),  one  obtains  the  following  expres¬ 
sions  among  the  TOEM: 

Cii2-f-4Cij5=|Cin,  Ci23-f  6Ci55-F8C4S«  =  0.  (10) 

In  order  to  examine  the  validity  of  the  “strong  Cauchy 
relations”  for  GaAs,  the  values  of  jCm,  C,,2-)-4C|j5,  and 
Ci23-(-6Ci44-j-8C«54  as  a  function  of  temperature,  evaluated 
from  the  data  shown  in  Fig.  4,  are  presented  in  Fig.  8.  As 
can  be  seen,  there  is  a  tendency  for  a  better  agreement  with 
the  Cauchy  relations  as  0  K  is  approached,  although  there 
appears  to  be  reasonable  agreement  at  room  temperature. 
Since  both  the  Keating  model  and  the  Cauchy  relations 
should  strictly  apply  only  at  0  K,  this  may  indicate  that  the 
lattice  interaction  in  GaAs  is  basically  of  the  central  type, 


100  200  irn 

TCMPERATURF  (K  ) 


FIG  4.  Third.order  coupling  consunt  A',  as  a  function  of  temperaiurr 


TABLE  VI.  Short-  and  long-range  inieraclion  contribution  to  the  roooi- 
lempcraturc  TOEM  (units  are  GPa) 


^lll 

C,,: 

^iii 

c,« 

^113 

Short  range 

-627 

-386 

-79 

42 

-219 

-41 

Long  range 

-I 

-  I 

-II 

-18 

-50 

-3 

and  the  deviations  from  the  Cauchy  relations  are  due  to 
thermal  effects.  An  analogous  observation  has  also  been 
made  for  germanium  and  silicon.' 


•4 

1 

< 


C.  Equation  of  state 

By  using  the  measured  values  of  the  TOEM,  an  equa¬ 
tion  of  state  (pressure- volume  relation)  for  the  material 
may  be  constructed.  One  of  the  simplest  and  most  straight¬ 
forward  is  the  Mumaghan  equation  of  stat^" 

y  !  p.-d/Bgl 

■  <"> 

where  Sq  is  the  bulk  modulus  and  Bq  its  pressure  deriva¬ 
tive,  both  at  P=0.  Bq  is  given  by 

Bi=-(l/3Bo)(C,-iC3).  (12) 

The  Mumaghan  equation  of  state  for  GaAs,  Eq.  (11),  is 
shown  in  Fig.  9.  The  dot  on  the  figure  is  the  experimentally 
determined  volume  at  the  onset  of  17.2  GPa  phase  transi¬ 
tion.  As  can  be  seen,  the  Mumaghan  equation  of  state 
predicts  a  somewhat  higher  compressibility  than  the  mea¬ 
sured  one. 


D.  Mode  GrUneisen  gamma  and  Grilnelsen  constant 

The  TOEM  are  closely  related  to  the  mode  Griineisen 
gammas  and  the  Griineisen  constant.’  Within  the  frame¬ 
work  of  the  anisotropic  continuum  model,’  the  mode  Grii¬ 
neisen  gammas  may  be  expressed  as 

r(f,N)  =  (Bo/2w)(pofV%  (13) 


II  MiT.KAll'RLiKi 

FIG  5  The  third-order  force  constanls  y.  6.  and  <  as  a  function  of 
temperature 
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FIG.  6.  The  three  TOEM,  Cm,  C|,j,  and  C,44  as  a  function  of  temper¬ 
ature. 


TEMPERATURE  (K) 


FIG.  8.  Measured  combinations  of  the  TOEM  compared  with  those  de¬ 
rived  from  the  Cauchy  relations,  as  a  function  of  temperature. 


Here,  N  is  a  unit  vector  in  the  propagation  direction  and  p 
the  polarization  index  (/>=  1,2,3).  For  a  cubic  material 
y(/>,N)  may  be  expressed  as 


r<;=  r(/>.N)C(p,N)(  1  jd£l  y(p.N)  j 


y(/’.N)  =  — (l/6io)(3£o+2t«-i-/:),  (14) 

where 

(15) 

fc(p.N)=c,i:,+Cii:2+(C,-C3)K3.  (i6) 

and 

Ki  (/>,N)  =N]U]+N\UI+N]UI 

K2ip,^)  =  (N2Ui+N3U2)^+(N3Ui+NiU2)^ 

+  iN,U2+N2UO\  (17) 

K2{p.N)=2{N2NjU2U3+N3NtU2U;+NtN2UiU2). 

The  mode  Griineisen  gammas  may  be  related  to  the  Grii- 
neisen  constant,^  Tc, 


In  this  expression,  C(j),fi)  is  the  heat  capacity  of  the  (p,N) 
mode,  n  the  spherical  angle,  and  the  integration  is  carried 
out  over  the  irreducible  part  of  the  Brillouin  zone,  i.e.,  the 
spherical  triangle  whose  apexes  are  [1(X)],  [1 10],  and  [1 1 1]. 

The  mode  ganunas  for  some  crystalline  directions  of 
high  symmetry  are  presented  in  Hg.  10.  It  is  noteworthy 
that  the  mode  gammas  for  the  slow  shear  mode  become 
negative  for  certain  directions.  This  raises  the  possibility 
that  the  thermal  expansion  will  become  negative  at  low 
temperatures,  since  this  is  the  lowcst-energy  mode,  and 
since  the  lowest-energy  modes  are  the  ones  excited  at  low 
temperatures.  This  is  borne  out  by  the  observation  that  the 
thermal  expansion  of  GaAs  becomes  negative  around 
40  K. 

In  the  limiting  cases  of  low  and  high  temperature  Eq. 
(18)  simplifies  significantly,  and  the  corresponding  values 
of  y  are  given  by 


TENfTER/MURE  'K/ 

FIG.  t.  I  be  thPBc  TOEKt,  Cjjj.  Cu^,  snd  «  function  of  temper- 

ilure. 


FIG.  9.  Munugban  equation  of  tute  for  GaAs.  The  date  point  is  an 
experimentally  determined  volume  at  the  onset  of  the  17.2  GPa  phase 
transition. 
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TABLE  VII.  Elastic  and  thermal  values  of  the  low.  and  high-temperature 
limits  of  the  Griineisen  constant. 
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FIG.  10.  GaAs  mode  Griineisen  gammas  for  some  crystalline  directions 
of  high  symmetry. 


Yl=  lS(p,N)-\ip,N){  2S(p,N)-^ 

AN  '  aN 

f  Ir(AN)](A^)-'.  (19) 

\  aN  ' 

where  the  summation  is  to  be  carried  out  over  all  modes,  N 
being  the  total  number  of  modes,  and  5(p,N)  the  sound 
velocity.  In  Table  VII,  the  values  of  and  as  deduced 
from  Eq.  ( 19)  are  compared  with  the  analogous  values  of 
the  thermal  Griineisen  constant,  evaluated  from  the  ther¬ 
mal  expansion  and  specific  heat.  As  can  be  seen,  there  is 
quite  good  agreement  between  the  elastic  and  thermal  y^, 
while  for  y£  the  elastic  value  is  somewhat  lower  than  the 
thermal  one.  But  in  this  case  the  comparison  is  not  very 
meaningful:  At  very  low  temperatures  the  thermal  y  rises 
steeply  as  the  temperature  is  lowered. 

V.  CONCLUSIONS 

The  present  work  demonstrates  that  the  Keating 
model  is  useful  in  the  evaluation  of  all  six  TOEM  of  QaAs 
at  room  temperature  and  shows  how  it  and  SHG  results 
can  be  used  to  evaluate  the  temperature  dependence  of  the 
TOEM.  The  agreement  between  experimentally  measured 
TOEM  at  room  temperature  and  the  ones  calculated  by 
combining  the  Keating  model  and  SHG  results  suggests 
that  the  Keating  model,  and  its  extension  to  include  long- 
range  Coulombic  interactions,  is  applicable  to  all  ciystals 
having  the  zinc-blende  structure.  We  suggest  that  the  same 
procedure  can  be  applied  to  other  zinc-blende  structure 


- — -C - 

Yl 

Yh 

Elastic 

0.38 

0.66 

Thermal 

0.60 

0.70  (  300  K) 

materials;  especially  the  other  III-V  compounds.  This  pro¬ 
duces  the  temperature  dependence  of  the  TOEM  of  III-V 
compounds  even  in  the  absence  of  direct  measurements  of 
all  of  them. 
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INTRODUCTION 

In  the  experiient  of  second  hanonic  generation 
(SHG),  the  capacitive  detector  is  used  to  Kasure  abso¬ 
lute  anplitude  of  acoustic  vave.‘^*  The  asseibly  of  the 
detector  and  saaple  can  be  siiplified  as  shovn  by  Pig.i. 
Usually  the  capacitive  detector  is  sechanical  dlsplace- 
nent  sensitive.  Vhen  the  saiple  is  piezoelectric,  hov- 
ever,  it  has  been  observed  that  capacitive  detector  gives 
output  even  there  is  no  DC-bias  applied  to  it. in  the 
present  paper,  the  origin  of  the  no  DC-bias  output  is 
discussed  and  its  influence  on  the  seasuresent  of  ultra¬ 
sonic  nonlinearity  paraneter  p  is  estisated.  The  calcula¬ 
tion  is  covared  with  eiperisent.  Although  the  analysis 
is  done  only  for  longitudinal  wave  along  Z-azis  of 
crystal  liRbO^  the  procedure  of  the  analysis  can  be 
easily  extended  to  any  piezoelectric  sediui  as  long  as 
there  exists  a  piezoelectric-stiffened  wve  in  certain 
direction. 


Moreover,  in  addition  to  incident  and  reflected  vaves 
there  should  be  a  reflected  evanescent  wve  A  in  the 
piezoelectric  sedius  in  order  to  fulfill  the  boundary 
conditions>^^  4^  satisfies  Eqs.U)  and  (2)  with  the 
soluticn  of  u^O,  and  4^(4. 

OPEK-CUCUIT 

In  open-circuit  case,  the  bottos  surface  of  sasple 
is  tnsetallized.  The  solution  for  reflected  evanescent 
wve  is 

M) 

here  b  is  penetration  depth  of  4. i.e.  4^-8  ihen  z^h.  The 
solution  of  one-diaensional  Laplacian  equation  is 

,{Z-d)  *4e 

vhicfa  also  satisfies  that  4:4^  vhen  z=d.  4^  is  the 
potential  developed  at  load  T.  use  of  boundary  conditions 
gives: 

(8.a) 

•I 

(8-C) 

In  addition,  the  circuit  equation 


ELECTRIC  FIELD  IN  AIRGAP 


Crystal  LiRbOjhas  sysKtry  of  3b.  The  longitudinal 
wve  along  Us  crystallographic  Z-axis  is  piezoelectric- 
stiffened.  Vithout  loss  of  generality  one  disensional 
problei  is  treated  here  because  nst  of  SHG  experisents 
are  perforsed  for  pure  longitudinal  vave  direction.  Hence 
the  coupling  equation  of  particle  displacesent  u  and 
electric  potential  4  can  be  expressed  as: 


should  be  applied.  Here,  A,  and  B,  are  avlitude  of 
incident  and  reflected  vaves,  respectively,  s  is  top 
surface  area  of  detector  button.  By  solving  Eqs(6)  the 
following  results  are  obtained; 


*•  l{ki-kdx^)**{kdx^)*]^'^ 


(7) 


Here  T^  and  0<  are  the  stress  and  electric  displacesent 
associated  with  nonlinear  constants  of  the  sedius. 

According  to  Eig.l.  there  is  sechanical  boundary 
condition  of  stress-free.  The  electric  boundsiy  condi¬ 
tions  consist  in  continuit,  nf  4  or  ronfr’clt:'  or  4 
D,f  depending  on  either  the  bottos  surface  of  the  saaple 
is  short  or  open  circuited.  Only  the  linear  part  of  the 
solution  is  given  here. 

In  Ue  airgap  between  8a^)le  and  detector  button 
there  exists  only  electric  potential  4^  which  satisfies 
Laplacian  Equation: 


•-a  -  *e  ^  •> ) 

*  T  'la  '  '  c  a 

l^^s  the  electroaechanical  coupling  factor. 


SBOBT-CIECUIT 


In  the  case  of  short-circuit,  the  bottos  surface  of 
saaple  is  coated  with  good  conductor  fila  and  grouided. 
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ABSTRACT 

Over  the  past  quarter-century  we  have  developed  techniques  for  evaluation  of  the 
nonlinear  properties  of  solids.  In  crystalline  solids  we  can  evaluate  third  order  elastic 
coi\stants  in  reliable  fashion.  Our  technique  makes  use  of  harmoitics  generated  by  the 
nonlinearity  of  the  meditun  as  finite  amplitude  ultrasonic  waves  propagate  Some  of  our 
recent  measurements  with  silicon  and  germanium,  two  diamond  lattice  solids,  are 
mentioned  by  way  of  illustration. 

Recently  we  studied  the  nonlinear  behavior  of  PZT  which  exhibits  a  very  large 
nonlinearity  parameter.  These  results  are  described  in  more  detail,  and  the  implication  of 
them  for  development  of  nonlinear  techniques  for  nondestructive  evaluation  of 
composites  is  stressed. 


^  I.  INTRODUCTION 

There  are  many  types  of  nonlinearity.  In  both  plasma  physics,  and  optics  one  can 
find  electromagnetic  nonlinearities.  The  nonlinearity  we  are  interested  in,  however,  is 
thermodynamical  or  mechanical  The  first  thermodynamical  (or  acoustical)  nonlinearity 
studied  probably  was  that  of  Poisson  in  1808.  Poisson^  studied  the  propagation  of  sound 
waves  of  finite  amplitude  in  an  ideal  gas.  The  first  study  of  the  propagation  of  ultrasonic 
waves  of  finite  amplitude  in  liquids  probably  was  that  of  Fox  and  Wallace^.  Subsequently 
Keck  and  Beyer  developed  a  theory  for  fluids^,  and  many  people  have  done  useful 
experiments^'  often  under  fimding  of  the  Navy.  The  foundation  of  the  study  of  the 
noiUinear  properties  of  solids  by  ultrasonic  techniques  was  laid  by  the  definition  of  third 
ord^  elastic  (TOE)  constants  in  the  1960's.  Since  that  time  in  my  laboratory  we  have  done  a 
number  of  experiments  that  suggest  to  us  that  noiUinear  properties  of  composites  might 
provide  the  basis  of  their  nondestructive  evaluation,  and  possibly  even  nondestructive 
evaluation  of  heat  damage  in  them. 

In  this  paper  I  propose  to  describe  pertinent  aspects  of  what  we  have  done  in  the 
past,  define  the  relationship  between  our  past  experience  and  characterization  of  the 
noiUinear  behavior  of  composites,  and  finally,  to  give  some  recent  resiUts  obtained  in  our 
laboratoiy  with  graphite  epoxy  composites. 

n.  THEORY 

A-  Linear  approximation-Hexagonal  Symmetry 

In  order  to  derive  the  linear  wave  equation  one  can  simply  define  the  elastic 
potential  energy  in  terms  of  the  elastic  constants: 

CijUTlijTlkl  .  0) 

ijkl 

This  strain  energy  substituted  into  Lagrange’s  equations  gives  the  linear  wave  equation  for 
principal  directions  a; 


j%  * 


where  K2  is  a  linear  combination  of  elastic  constants.  Our  experience  leads  us  to  assume 
that  epoxy  composites  will  be  describable  in  terms  of  equations  appropriate  to  hexagonal 
symmetry.  To  evaluate  the  elastic  constants  of  composites  we  propagate  ultrasonic  waves 
along  the  X3  axis.  In  this  case 

pv^fOOl)  =  C33  for  longitudinal  waves  G) 

and 

pv^(lOO)  =  C44  for  transverse  waves.  (4) 

For  propagation  along  any  direction  perpendicular  to  the  X3  axis 

pv^ClOO)  =  Cii  for  longitudinal  waves  (5) 

and 

pv^(lOO)  «  C44  for  tnmsverse  waves  polarized  in  the  (001)  direction  (6) 

pv^dOO)  =  2  ^^11 "  ^12)  for  transverse  waves  polarized  in  the  (010)  direction  (7) 

6.  Nonlinear  approximation  -  Cubic  Symmetry 

Most  of  our  experiments  on  the  propagation  of  finite  amplitude  ultrasoiuc  waves 
have  been  done  with  solids  of  cubic  symmetry,  although  we  have  worked  out  the  nonlinear 
theory  for  any  syiiunetry^-  To  describe  our  experiments,  then,  it  is  necessary  only  to  state 
theoretical  backgroimd  for  cubic  symmetry.  We  begin  again  with  the  elastic  potential 
energy,  but  now  keep  higher  order  terms: 

4Kn)  ~  ^  Cijiciinn  TlijTjuTJjjuj  +  •  •  •  (8) 

ijkl  ■  ijklmn 

Now  we  define  the  Lagrangian  function 


(9) 


and  substitute  into  Lagrange's  equations: 


(10) 


The  result  is  that  for  prindpcil  directions  in  a  cubic  crystal  we  can  write  the  nonlinear  wave 
equation  in  the  form 


Po 


+  (3K2  +  K3) 


au 

aa  aa2 


(11) 


where  expressions  for  K2  and  K3  for  a  cubic  lattice  are  given  in  Table  I. 


Table  I.  K2  and  K3  for  principal  directions  in  a  cubic  crystal 


Direction  K2  K3 


1100] 

Cll 

an 

[110] 

Cii  +  C12  +  2C44 

Cm  +  3Cii2  +  12Ci66 

2 

4 

The  solution  of  Eq.  (11)  pertinent  to  the  present  discussion  is  obtained  under  the 
assumption  that  a  sinusoidal  disturbance  is  generated  at  a  =  o.  At  a  distance  a  from  the 
sinusoidal  driver  the  solution  takes  the  form 


U  =  Ai  sin(ka  -  cot) 


3K2  + 


K2 


Aj  k^a  cos  2(ka  -  cot)  +  • 


(12) 


Our  procedure  involves  the  measurement  of  the  amplitude  of  the  fundamental  Ai  and  the 
amplitude  of  the  second  harmonic 


A2  =  -^l^^A2k2a 
K2 

The  nonlinearity  parameter  is  defined  as 


(13) 


p_  3K2  +  K3 
K2 


(14) 


From  Eqs.  (13)  and  (14)  we  can  define  the  nonlinearity  parameter  in  terms  of  measured 
quantities  as 


In  order  to  determine  the  nonlinearity  parameter,  then,  we  need  only  to  measure  the 

amplitudes  of  the  fundamental  and  that  of  the  second  harmonic  A2.  The  propagation 

constant  k  =  ^  can  be  determined  from  the  frequency  and  the  sound  velocity.  The 

X 

quantity  a  is  the  sample  length. 


HL  MEASUREMENT  TECHNIQUE 

Measurement  of  the  amplitudes  of  the  fundamental  and  the  second  harmonic  can  be 
made  with  a  variation  of  the  capacitive  microphone  as  shown  in  Fig.  1.  The  detector  button 
produces  an  air  gap  of  the  order  of  lOp  between  the  button  and  the  end  of  the  sample.  This 
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Figure  1.  Capacitive  detector  for  measurement  at  room  temperature. 

parallel  plate  capacitor  is  very  sensitive.  It  allows  us  to  measure  amplitudes  as  small  as  10*^ 
Angstroms.  The  detector  is  used  in  the  apparatus  shown  in  Fig.  2  for  the  measurement  of 
fundamental  and  second  harmonic  amplitudes. 


IV.  RESULTS 

Measxired  nonlinearity  parameters  as  a  function  of  temperature  are  shown  in  Fig.  3. 
The  temperature  dependence  is  relatively  small  so  one  can  correlate  different  crystalline 
structures  and  bonding  with  the  magnitudes  of  the  nonlinearity  parameters.  In  most 
materials  the  magnitudes  of  the  , nonlinearity  parameters  are  between  2  and  15.  Table  n 


Figure  2,  Experimental  arrangement  for  measurement  of  the  nonlinearity  parameter. 
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Figure  3.  Measured  temp>erature  dependence  of  nonlinearity  parameters. 


Table  IL  Compauison  of  ultrasonic  nonlinearity  parameters 


MA'^'ERIAL 

OR  STRUCTURE 

BONDING 

Pav 

22ncblend 

Covalent 

2.2 

Fluorite 

Ionic 

3.8 

PCC 

Metallic 

5.6 

KX  (Inert  gas) 

Van  der  Waals 

6.4 

BCX 

Metallic 

8.2 

NaQ 

Ionic 

14.6 

Fus^  Silica 

Isotropic 

-3.4 

YBa2Cu307.6 

Ceramic 

Isotropic 

143 

gives  a  correlation  between  average  nonlinearity  parameters  and  crystalline  structure. 
Zincblend,  or  diamond  lattice  solids  with  covalent  bonding,  have  the  smallest  nonlinearity 
parameters  of  about  2,  and  Nad  with  ionic  bonding  has  the  largest  noidinearity  parameter, 
14.6.  The  difference  between  nonlinearity  parameters  is  real.  Therefore  variations  of 
nonlinearity  parameter  may  be  correlated  with  material  behavior  and  ultimately  may  serve 
as  basis  of  a  nondestructive  evaluation  technique  for  these  materials,  and  possibly  others 
such  as  composites. 


A.  Correlation  of  results  for  diamond  lattice  solids 

In  order  to  test  specific  results  on  diamond  lattice  solids  we  evalxiated  nonlinearity 
parameters  of  silicon  and  germanium.  From  Table  I  we  were  able  to  isolate  specific 
combinations  of  third  order  elastic  constants.  The  combinations  we  could  isolate  are  shown 
in  Figs.  4  and  5.  The  fact  that  the  temperature  behavior  of  these  two  diamond  lattice  solids 
is  similar  indicates  that  our  data  are  meaningful.  In  addition,  we  were  able  to  evaluate  the 
behavior  of  the  third  order  elaistic  constants  of  a  cential  forces,  nearest-neighbor  model.  For 
such  crystals  C112  +  4  C166  =  5/2  Cm  and  C123  =  C 144  =  C456  =  0.  Examination  of  Figs.  4  and  5 
reveals  that  such  relationships  are  approached  at  absolute  zero  of  temperature  for  both 


TEMPERATUR£(*K) 


to 

8  / 

uj-20 

O 

t- 


Rguie4.  Thizd  order  elastic  comtants 
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Hgure  5.  Third  order  elastic  constants 
of  gennaidum. 


silicon  and  germanium  .  This  is  consistent  with  independent  information  about  the  two 
crystals. 

A  further  test  of  the  validity  of  the  data  can  be  made  if  all  six  third  order  elastic 
constants  can  be  isolated  from  the  data.  We  isolated  all  six  third  order  elastic  constants, 
then  used  them  to  calculate  the  Griineisen  parameter  which  can  be  determined 
independently  from  thermal  expansion  data.  The  comparisor^  are  shown  in  Figs.  6  zmd  7. 
The  minimum  that  shows  up  in  the  thermal  experiment  at  approxinuitely  0.1  6  (where’ 6  is 
the  Debye  temperature)  is  reproduced  by  our  data  at  the  proper  temperature.  Note  that  our 
data  are  better  than  those  of  other  experimenters  in  this  respect.  One  can  conclude, 
therefore,  that  the  magnitudes  of  the  nonlinearity  parameters  can  be  correlated  with  real 
physical  behavior  of  silicon  and  germanium. 
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Figure  6.  Griineisen  parameter  of  silicoa 


Rgure  7.  .Grimeisen  parameter  of 
germanium. 


B.  Effect  of  Material  Behavior  on  Nonlinearity  Parameter 

The  effect  of  material  behavior  on  measured  values  of  the  nonlinearity  parameter  can 
be  illustrated  by  recent  measurements  in  PZT.  As  is  well-known  neeu"  the  Curie 
temperature  PZT  goes  from  the  room-temperature  symmetry  (which  can  be  either 
rhombohedial  or  tetragonal,  dej>ending  on  the  ratio  of  Pb'n03  to  PbZr03)  to  cubic  symmetry. 
This  atomic  rearrangement  is  somewhat  arudogous  to  what  happens  in  metals  during 
work-harderring. 

Measured  values  of  the  nonlinearity  parameters  of  two  samples  of  PZT  are  shown  in 
Figures  8  and  9.  In  Fig.  8  the  Si  sample  goes  from  rhombohedral  to  cubic  at  the  Curie 
temperature.  In  this  case  the  nonlinearity  parameter  goes  from  p  •=  4  at  room  temperature 
to  p  *  250  at  the  Curie  temperature.  In  Fig.  9  the  nonlinearity  parameter  of  K1  PZT  is  shown 
to  go  from  its  room  temperature  value  of  p  *  8  to  p  *  1500  at  the  Curie  temperature  as  the 
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sample  goes  front,  tetragonal  to  cubic  symmetry.  Such  large  variatioi\s  in  the  nonlinearity 
parameter  are  associated  with  molecular  rearrangements  in  the  materials  themselves.  Such 
rearrangements  also  may  take  place  in  Heat  damage  to  composites.  This  hypothesis 
suggests  that  measurement  of  the  nonlinearity  parameter  of  composites  might  be  a  very 
Useful  experimental  program.  It  could  result  in  a  completely  new  nondestructive 
evaluation  technique.  Such  a  technique  would  depend  on  variation  in  nonlinear 
properties  of  composites,  information  that  heretofore  has  been  ignored. 

C  Experimental  results  with  composites. 

As  a  preliminary  set  of  experiments  we  have  measured  the  velocity  of  ultrasonic 
waves  in  graphite  epoxy  rods.  The  results  are  shown  in  Table  EL  It  is  to  be  noted  that  the 
velocity  of  ultrasonic  waves  along  the  fiber  direction  (axial  velocity)  is  three  times  the 
velocity  along  a  perpendicular  direction  (radial  velocity).  If  one  assumes  that  the  hexagonal 
model  is  appropriate,  then  this  velocity  difference  translates  into  a  factor  of  10  in  the  ratio  of 
the  elastic  cortstants,  since 

Cij  =  pV2  . 

Such  velocities  should  be  measured  with  samples  that  actually  are  subjected  to  heat 
damage.  A  first  test  of  the  value  of  the  hexagonal  model  of  graphite  composites  would  be  to 
determine  the  result  of  heat  damage  on  this  ratio  of  velocities. 


Table  EL  Measured  ultrasonic  velocities  and  elastic  constants  of  graphite-epoxy  composites. 


PROPAGATION 

COMPRESSIONAL 

■  ■  "  ■  "  ■ 

ELASTIC 

DEJECTION 

VELOOTY 

CONSTANT 

(M/S) 

(DYN/Qd2) 

AXIAL 

8088  89 

C33  =  (1.074  ±  0.031)  xl0l2 

RADIAL 

2S35±59 

Cii  =  (0.1046  ±  0.0056)  x  lO^^ 

p  *  1.627  ±0.011  g  /  cm^ 


D.  Nonlinear  measurements  in  composites 

Although  experimental  research  is  yet  to  be  done  on  the  nonlinear  properties  of 
composites/  a  theoretical  basis  for  it  exists.  We  have  calculated  the  behavior  of 
contributions  to  the  nonlinear  behavior  of  crystals  of  all  symmetries.  For  hexagoxuil 
symmetry  we  find  that  the  velocity  is  the  same  in  all  direction^  in  the  basal  plane.  In 
contrast/  our  calculations  show  the  behavior  of  the  third  order  elastic  constants  in  the  basal 
plane  given  in  Fig.  10.  For  the  [100]  direction  the  appropriate  third  order  elastic  constant  is 
Cm.  For  the  [010]  direction  it  is  C222*  Since  in  general  C222  ^  Qll  for  hexagonal  symmetry/ 
the  six-fold  symmetry  is  in  evidence  in  the  nonlinear  behavior  of  a  hexagonal  crystal  even 
though  it  may  be  imobservable  in  the  velocity/  a  linear  phenomenon. 


C010] 


'v.  CONCLUSION 


Variations  in  the  nonlinearity  parameter  already  have  been  correlated  with  ultimate 
yield  strength  of  cubic  solids,  with  hardness  in  steels,  with  crystalline  structure  and  bonding, 
and  with  thermal  expansion  coefficient  Some  experimental  results  taken  in  our  laboratory 
are  discussed  and  are  shown  to  indicate  that  measurement  of  the  nonlinear  properties  of 
composites  could  result  in  a  very  sensitive  technique  for  nondestructive  evaluation  of 
composites.  Such  a  technique  might  be  sensitive  enough  to  detect  heat  damage  in 
composites  lor\g  before  it  can  be  detected  by  other  tedmiques. 
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